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Abstract 



We provide a complete and rigorous description of phase transitions for ki- 
netic models of self-propelled particles interacting through alignment. These 
models exhibit a competition between alignment and noise. Both the align- 
ment frequency and noise intensity depend on a measure of the local align- 
ment. We show that, in the spatially homogeneous case, the phase transition 
C^ . features (number and nature of equilibria, stability, convergence rate, phase 

diagram, hysteresis) are totally encoded in how the ratio between the align- 
ment and noise intensities depend on the local alignment. In the spatially 
inhomogeneous case, we derive the macroscopic models associated to the sta- 
ble equilibria and classify their hyperbolicity according to the same function. 
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1 Introduction 

In this work we provide a complete and rigorous description of phase transitions 
in a general class of kinetic models describing self-propelled particles interacting 
through alignment. These models have broad applications in physics, biology and 
social sciences for instance for the description of animal swarming behavior or opin- 
ion consensus formation. Their essential feature is the competition between the 
alignment process which provides self-organization, and noise which destroys it. An 
important point is that both the alignment frequency and noise intensity depend 
on a measure of the local alignment denoted by \J\. The phase transition behavior 
in the spatially homogeneous case is totally encoded in the ratio between these two 
functions denoted by k{\J\). Namely we have the following features: 

(i) The function k gives rise to an algebraic compatibility relation whose roots 
provide the different branches of equilibria of the kinetic model. One distin- 
guished branch is given by isotropic or uniform distributions which correspond 
to no alignment at all, i.e. \J\ = 0. The other branches are associated to 
non-isotropic von Mises-Fisher distributions associated to non-zero \J\. 

(ii) The stability of these various equilibria is completely determined by the mono- 
tonicity of a function derived from k around these roots and there exists an 
exponential rate of local convergence of the solution to one of these stable 
equihbria. 

(iii) The global shape of this function k provides the phase diagram which encodes 
the order of the associated phase transitions. According to its monotonicity, 
these can be second-order phase transitions, first-order phase transitions with 
hysteresis behavior or even be more complex. For second-order phase tran- 
sition, we give an explicit formula for the critical exponent in terms of the 
local behavior of k. The involved phase transitions are spontaneous symme- 
try breaking phase transitions between isotropic and non-isotropic equilibria. 
Such phase transitions appear in many branches of physics, such as spon- 
taneous magnetization in ferromagnetism, nematic phase transition in liquid 
crystals and polymers, Higgs mechanism of mass generation for the elementary 
particles. 



(iv) In the spatially inhomogeneous case, we can derive the hydro dynamic equa- 
tions associated to both the isotropic and non-isotropic stable equilibria (the 
former leading to diffusion behavior, the latter to hyperbolic models). The 
hyperbolicity is again completely determined by this function, and is linked to 
the critical exponent in the case of a second-order phase transition. 

To our knowledge, this is the first time that a complete mathematical theory of 
phase transitions in a physics system can be rigorously derived and related to one 
single object with high physical significance: this function k. One of the main 
achievement of this work is Theorem [2|, which provides part of point (ii) above, 
namely the nonlinear stability of the non-isotropic equilibria (the von Mises-Fisher 
distributions) when the function associated to k is increasing. To be more precise, 
let us write this set of equilibria as {Jq, ^2 G S} (it has the same symmetries as 
the unit sphere S of M"", n being the dimension of the model). Then, we have a 
rate of convergence A and two positive constants S and C such that, if the initial 
condition /o satisfies ||/o — f^W < S for some fi G §, then there exist fioo G § such 
that for all t > 0, the solution / of the spatially homogeneous model satisfies 

ll/W-/^LIKC'||/o-/^^||e-^*. 

This stability result takes place in the Sobolev space H^ as long as s > ^^. In 
previous works (in the case where the function k is linear) such as [T7j or [IB] (for 
the Kuramoto model in dimension n = 2, where a precise study of the attractor is 
performed), the exponential convergence with rate /3 was only proven for all /3 < A, 
and the existence of such a constant C independent of /o was lacking. 

Self-propelled particle systems interacting through alignment have been widely 
used in the modeling of animal swarms (see e.g. the review [21] and [2| [^[5| [271125] ) . 
Kinetic models of self-propelled particles have been introduced and studied in [3],|ll 
[T21l2U[|21j . Here, specifically, we are interested in understanding phase transitions and 
continuum models associated to the Vicsek particle system [2H]- Phase transitions 
in the Vicsek system have been widely studied in the physics literature [HE] . There 
have been some controversies whether the involved phase transitions were first or 
second order. In some sense, this paper provides a complete answer to this question, 
at least in the kinetic framework. 

The passage from the kinetic to macroscopic descriptions of the Vicsek system 
has first been proposed in [12]. Further elaboration of the model can be found 
in [TTlfTB]. The resulting continuum model is now referred to as the Self-Organized 
Hydrodynamic (SOH) model. In these derivations of the SOH, the noise and align- 
ment intensities are functions of the local densities and not of the local alignment. 
No phase transition results from this choice but the resulting SOH models are hy- 
perbolic. In [iniE], alignment intensity has been made proportional to the local 
alignment. Second-order phase transition have been obtained. However, the re- 
sulting SOH model is not hyperbolic. In the present paper, we investigate general 
relations between the noise and alignment intensities upon the local alignment \J\. 
As described above, the phase diagram becomes extremely complex and its com- 
plexity is fully deciphered here. The kind of alignment phase transition that we find 
here is similar to nematic phase transitions in liquid crystals, polymer dynamics and 
ferromagnetism [7l|Tll[I51[23l[21] • 



The organization of the paper is as follows. In section El we derive the kinetic 
model from the particle system and determine its equilibria. In section E], we study 
the stability of these equilibria in the spatially homogeneous case and find the rates 
of convergences of the solution to the stable ones. Then, in section HJ we use these 
results to study two examples respectively leading to second order and first order 
phase transitions, and in the case of first order phase transitions, to the hysteresis 
phenomenon. Finally, in section we return to the spatially inhomogeneous case 
and investigate the macroscopic limit of the kinetic model towards hydrodynamic 
or diffusion models according to the considered type of equilibrium. For the hy- 
drodynamic limit, we provide conditions for the model to by hyperbolic. Finally, a 
conclusion is drawn in section [HI We supplement this paper with appendix |X] which 
provides elements on the numerical simulation of the hysteresis phenomenon. 

2 Kinetic model and equilibria 

In this section, we derive the mean-field kinetic model from the particle system, and 
determine its equilibria. We begin with the particle model in the next section. Then, 
in section [22] we derive the mean-field limit. The space-homogeneous case will be 
highlighted in section 12.31 and the equilibria will be determined in section 12.41 

2.1 The particle model 

We consider a system of a large number N of socially interacting agents defined 
by their positions Xi G M" and the directions of their velocities cuj G S (where S is 
the unit sphere of M*^). They obey the following rules, which are a time continuous 
version of those of the Vicsek model [2H] : 

- they move at constant speed a, 

- they align with the average direction of their neighbors, as a consequence of 
the social interaction. 

- the directions of their velocities are subject to independent random noises, 
which expresses either some inaccuracy in the computation of the social force 
by the subject, or some trend to move away from the group in order to explore 
the surrounding environment. 

These rules are expressed by the following system of stochastic differential equations: 

dXi = aUidt, (2.1) 



do;, = z/(|j;|)P^x iJidt + J2t{\J,\)P^^ o dBl (2.2) 

7- a ^ 

^^ = rw:v J^ = ^Y.^<{\^^-U)^i■ (2.3) 

Eq. (12. Ij) simply translates that particle i moves with velocity awj. The first term at 
the right-hand side of 02.21) is the social force, which takes the form of a relaxation of 



the velocity direction towards the mean direction of the neighbors Ui, with relaxation 
rate u (the operator P ± is the projection on the tangent space orthogonal to Ui, 
ensuring that Ui remains a unit vector). Eq. (12.31) states that the mean direction is 
obtained through the normalization of the average current J7i, itself computed as the 
average of the velocities of the particles. This average is weighted by the observation 
kernel K, which is a function of the distance between the test particle i and its 
considered partner i. Without loss of generality, we can assume that / -ft'(lCl) d,^ = 1. 
The second term of (12.21) models the noise in the velocity direction. Eq. (12.21) must 
be understood in the Stratonovich sense (as indicated by the symbol o), with A^ 
independent standard Brownian motions Bl in M". The quantity r > is the 
variance of the Brownian processes. 

In this paper, we assume that the relaxation rate u and the noise intensity r 
are functions of the norm of the current \J'\. The present hypothesis constitutes a 
major difference with previous works. Indeed, the case where u and r are constant 
has been investigated in [12], while the case where z^(|i7|) = \J'\ and r = 1 has been 
treated in [TU]. We recall that no phase transition appears at the macroscopic level 
in the first case while in the second case, a phase transition appears. This phase 
transition corresponds to a change in the number of equilibria as the density crosses 
a certain threshold called critical density. The critical exponent is 1/2 in this case. 
Here, we investigate the more general case of almost arbitrary dependences of u 
and r upon \J'\, and show that the phase transition patterns can be much more 
complex than those found in [TU]. For later convenience, we will denote by Tq > 
the value of r(0). 

To understand why \J'\ is the crucial parameter in this discussion, let us introduce 
the local density pi and order parameter (or mean alignment) q as follows: 

17-1 1 ^ 



a Pi iV ^ 

and we note that ^ q ^ 1. The value Cj ~ corresponds to disorganized motion, 
with an almost isotropic distribution of velocity directions, while Q ~ 1 characterizes 
a fuUy organized system where particles are all aligned. Therefore | J7i| appears as the 
"density of alignment" and increases with both particle density and order parameter. 
This paper highlights that the dependence of u and r upon \J'i\ acts as a positive 
feedback which triggers the phase transition. Besides, in [TH], it has been shown that 
making u and r depend on the density p only does not produce any phase transition, 
and that the recovered situation is qualitatively similar to that of [12] • The present 
work could be extended to u and r depending on both p and \J'\ at the expense 
of an increased technicality, which will be omitted here. The present framework is 
sufficient to cover all interesting situations that can be desirable at the macroscopic 
scale. 

2.2 Mean-field derivation of the kinetic model 

The first step in the study of the macroscopic behaviour of this system consists 
in considering a large number of particles. In this limit, we aim at describing the 
evolution of the density probability function f^{x,u},t) of finding a particle with 



direction u at position x. This has been studied in [3] in the case where i^dJ/"!) = \J'\ 
and r = 1. It is nearly straightforward to perform the same study in our more general 
case. 

For convenience, we will use the following notation for the first moment of a 
function / with respect to the variable u (the measure on the sphere is the uniform 
measure such that /§di^ = 1): 



Jf{x,t) = Lo f{x,uj,t)duj. (2.4) 

For the following, we will assume that: 
Hypothesis 2.1. 

(i) The function K is a Lipschitz bounded function with finite second moment. 

(a) The functions \J\ i— ?■ |j| and \J\ H- t(| J|) are Lipschitz and bounded. 

In these conditions the mean-field limit of the particle model is the following 
kinetic equation, called Kolmogorov-Fokker-Planck equation: 

dj + au- VJ + iy{\Jf\)V^ ■ (P^^uf f) = r{\Jf\)A^f (2.5) 

with 

Jf{x,t) = a{K*Jf){x,t), ^f = ^y (2.6) 

where * denotes the convolution in M" (only on the x variable), A^ and V^j- stand 
for the Laplace-Beltrami and divergence operators on the sphere §. 
More precisely, the following statements hold: 

Proposition 2.1. ///q is a probability measure on M" x§ with finite second moment 
in X G M", and if (X°, w°)jg|i7v| a'^e A^ independent variables with law /o, then: 

(i) There exists a pathwise unique global solution f to the particle system (12. ip - 
(12. 3p with initial data (X°,u;°). 

(a) There exists a unique global weak solution of the kinetic equation (12.50 with 
initial data /q. 

(Hi) The law f^ at time t of any of one of the processes {Xi,u}i) converges to f 

as N ^ oo. 

The proof of this proposition follows exactly the study performed in [5], using 
auxiliary coupling processes as in the classical Sznitman's theory (see [25]), and is 
omitted here. Let us make some comment on the structure of the kinetic equa- 
tion (12. 5p . The first two terms of the left hand side of (12.50 correspond to the free 
transport with speed given by au. It corresponds to (12.10 in the particle model. The 
last term of the left hand side corresponds to the alignment mechanism towards the 
target orientation Uf, with intensity z^(|j7/|), while the term at the right hand side 
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is a diffusion term in the velocity variable, with intensity t(|J7/|). These two terms 
correspond to 02.21) in the particle model. We will see in (12. 7p and f l5.6p that these 
two terms, under certain assumptions (spatially homogeneous case, or expansion in 
terms of a scaling parameter rj), behave as a local collision operator Q, only acting 
on the velocity variable u. Finally, the convolution with K in (12.61) expresses the 
fact that i7/ is a spatial averaging of the local momentum Jf defined in (12. 4p . it 
corresponds to the definition (12.30 in the particle model. 



2.3 The space-homogeneous kinetic model 

The hydrodynamic limit involves an expansion of the solution around a local equi- 
librium (see section EH]). Therefore, local equilibria of the collision operator Q are 
of key importance. We will see that such equilibria are not unique. The existence 
of multiple equilibria requires an a priori selection of those equilibria which make 
sense for the hydrodynamic limit. Obviously, unstable equilibria have to be ignored 
because no actual solution will be close to them. In order to make this selection, in 
the present section, we consider the spatially homogeneous problem. To the most 
possible exhaustive way, in section El we will determine the stable equilibria and 
characterize the convergence rate of the solution of the space-homogeneous problem 
to one of these equilibria. In section HJ we will illustrate these results on two exam- 
ples. Finally, in sectional we will deal with the spatially non- homogeneous case and 
apply the conclusions of the spatially homogeneous study. 

The spatially homogeneous version of this model consists in looking for solu- 
tions of the kinetic equation (12. 5p depending only on u and t. Obviously, such 
solutions cannot be probability measures on M" x § any more, so we are looking 
for solutions which are positive measures on §. In that case, J7/ = aJf, and (up to 
writing i^(| J/|) = i^(a| J/|) and t(| J/|) = T(a| J/|)) the kinetic equation (12.50 reduces 
to 

dtf = Q{f), (2.7) 

where the operator Q is defined by 

Qif) = -^(l^/l)V. ■ (P^xf]; /) + ri\Jf\)A^f, (2.8) 

where Qj = -ppT and where we have dropped the "hats" for the sake of clarity. Let 
us remark that by hypothesis 12. 1^ we do not have any problem of singularity of Q 
as \Jf\ — )• 0: if | J/| = 0, we simply have Q{f) = TqA^/. 

The investigation of the properties of the operator Q is of primary importance, 
as we will see later on. For later usage, we define 

H\J\) = '^y Hr)=rk{s)ds, (2.9) 

r(|J|) Jo 

so that $(|J|) is an antiderivative of k: 4% = /i;(|J|). The space-homogeneous 
dynamics corresponds to the gradient flow of the following free energy functional: 

T{f)= f f\nfdio-mJf\). (2.10) 

Js 



Indeed, if we define the dissipation term X'(/) by 

W) = r{\Jf\) f f |V.(ln/ - k{\Jf\) u ■ fi^)pda;, (2.11) 

we get the following conservation relation: 

^^(/) = -P(/) ^ 0. (2.12) 

The main ingredient to derive this relation is the identity P^i^Vtf = V^(cj ■ f^/)- 
Therefore, the collision operator Q defined in 02.81) can be written: 

Qif) = ri\Jf\)V^ ■ [/ V.(ln/ - A;(|J^|)c^ ■ Qf)]. (2.13) 

Finally, since 

^T= fdtf{\nf-k{\Jf\)u:-nf)dio, 
at Js 

using (12.71) . (I2.13P and integrating by parts, we get (I2.12p . 

We first state results about existence, uniqueness, positivity and regularity of 

the solutions of (12. 7|) . Under hypothesis 12.11 we have the following 

Theorem 1. Given an initial finite nonnegative measure /o in H''^{S), there exists a 
unique weak solution f of (12.71) such that /(O) = /q. This solution is global in time. 
Moreover, f e C\R%,C°^{S)), with f{uj,t) > for all positive t. 

Finally, we have the following instantaneous regularity and uniform houndedness 
estimates (form G N, the constant C being independent of fo): 

\\m\\u^^c(i + j^)\\fo\\is. 

The proof of this theorem follows exactly the lines of the proof given in [T7j 
for the case where z^(|J|) = |J|, and will be omitted here. Let us remark that 
here we do not need the bounds on |j| and on r provided by hypothesis 12.11 
since the positivity ensures that \J\ takes values in [0,po]) where po is the total 
mass of /o (a conserved quantity). Therefore r is uniformly bounded from below in 
time, by a positive quantity Tmm, and |j| is also uniformly bounded from above 
in time. Finally, the fact that / is only C^ in time comes from the fact that the 
proof only gives / G C([0,T], if*(§)) for all s, and we use the equation to get one 
more derivative. We could obtain a better time regularity at the price of a better 
regularity for the functions |j| and on r. 

2.4 Equilibria 

We now define the von Mises-Fisher distribution which provides the general shape 
of the non-isotropic equilibria of Q. 

Definition 2.1. The von Mises-Fisher distribution of orientation f2 G S and con- 
centration parameter k, ^ is given by: 



The order parameter c{k) is defined by the relation 

JM^n = c(«)^, (2.15) 

and has expression: 



cw = ^°::°,::.„...:r, "" . (2.16) 



The function c : k G [0, oo) H- c{k) G [0, 1) defines an increasing one-to-one 
correspondence. Tlie case k = c(k) = corresponds to the uniform distribution, 
while when k is large (or c{k) is close to 1), the von Mises-Fisher distribution is 
close to a Dirac delta mass at the point ^2. 

For the sake of simplicity, we will assume the following: 

Hypothesis 2.2. The function \J\ i— )• A;(|J|) = ^1L| is an increasing function. We 
denote by j its inverse, i.e. 

K = k{\J\)^\J\=3{K). (2.17) 

This assumption is not critical. It would be easy to remove it at the price of 
an increased technicality. Additionally, it means that when the alignment of the 
particles is increased, the relative intensity of the social force compared to the noise 
is increased as well. This can be biologically motivated by the existence of some 
social reinforcement mechanism. It bears analogies with Diffusion Limited Aggrega- 
tion (see [30]), in which the noise intensity is decreased with larger particle density. 
This can also be related with what is called "extrinsic noise" in [T], where the noise 
corresponds to some uncertainty in the particle-particle communication mechanism. 
Indeed in this case, the intensity of the noise increases when \J\ decreases. 

The equilibria are given by the following proposition: 

Proposition 2.2. The following statements are equivalent: 

(i) f G C^iS) and Qif) = 0. 

(ii) f eC\S) andVif) =0. 

(Hi) There exists p ^ and fi G S such that f = pM^n, where k ^ satisfies the 
compatibility equation: 

j{k) = pc{k). (2.18) 

Sketch of the proof. The proof is identical to that of [17], and we just summarize 
the main ideas here. The main ingredient is to observe that Q{f) (or T>{f)) is equal 
to zero if and only if / is proportional to Mk(\jj:\)nf This is quite straightforward 
for V using (12. 111) . For Q, it follows from the following expression: 



W) = ^(l^/I)v. 



'^■"""W. 



HUfD^f 



:2.i9) 



This expression comes from Definition 12.11 which gives first 

1 \ -MI>//l)V.(a;-fi/) k{\Jf\) 



Mki\jf\)nfJ Mk{\j^.\)^j Mk{\jj\)^j 



and therefore, applying the chain rule to the right-hand side of (12.191) . we recover 
the definition of Q given in f l2.8p . Hence, we obtain 



JQU) 



f 



MkiijfDQf 



duj = -r{\Jf\) f 



f 



Mk(\jf\)nf 



M^j^l^Ufduj. 



So if Q{f) = 0, we get that jj — is equal to a constant. Conversely if / is 

proportional to M^dj^Dn^, we directly get with (12.191) that Q{f) = 0. 

Now if / is proportional to M^dj^Dn^, we write / = pMf^Q, with n = k{\Jf\), 
which corresponds to |J/| = j{k) thanks to (12.171) . But then by (I2.15p . we get 
that I J/l = pc{K,), which gives the compatibility equation (12.181) . Conversely, if we 
have {iii), we also get that \Jf\ = pc{k) = j{k) and so k = k{\Jf\), which gives 
that / is proportional to Mfc(|j^.|)Q^. D 

We now make comments on the solutions of the compatibility equation (12.181) . 
Let us first remark that the uniform distribution, corresponding to k = is always 
an equilibrium. Indeed, we have c(0) = j(0) = and (12.181) is satisfied. However, 
Proposition 12.21 does not provide any information about the number of the non- 
isotropic equihbria. The next proposition indicates that two values, p* and pc, that 
can be expressed through the function k only, are important threshold values for the 
parameter p, regarding this number of non-isotropic equilibria. 

Proposition 2.3. Let p > 0. We define 

p, = hm 44 = lim „'f'„, = lim -^, (2.20) 

^ «^o c{k) \J\^o c{k{\J\)) |J|->oA;(|J|)' ^ ^ 

P* = inf 44 = inf ^77^, (2.21) 

^e{o,K^a^) c{k) \J\>o c{k{\J\)) ^ ' 

where pc > may be equal to +oo, where K^ax = ^™^\j\^ook{\J\) , and where we 
recall that n denotes the dimension. Then we have p^ ^ p^, and 

(i) If p < p^, the only solution to the compatibility equation is n = and the only 
equilibrium with total mass p is the uniform distribution f = p. 

(a) If p > p^, there exists at least one positive solution k, > to the compatibility 
equation (12.181) . It corresponds to a family {pMi^Q,Q G §} of non-isotropic 
equilibria. 

(iii) The number of families of nonisotropic equilibria changes as p crosses the 
threshold pc (under regularity and non- degeneracy hypotheses that will be pre- 
cised in the proof, in a neighborhood of pc, this number is even when p < pc 
and odd when p > pc). 

Proof. Some comments are necessary about the definitions of pc and p*. First note 
that, under hypotheses 12.11 and 12.21 k is defined from [0,+oo), with values in an 
interval [0,Kmax), where we may have Umax = +cxd. So j is an increasing func- 
tion from [0,Krnax) outo IR+, and this gives the equivalence between the two terms 
of (12.211) . Thanks to hypothesis 12. 1^ we have 

A;(|J|) = ^|J|+o(|J|)as |J|^0, 

To 
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with To = r(0) and Ui = lim|j|_»o ^\j, , and the last term of f l2.20p is well defined 
in (0, +cxd] (we have pc = ^ ii i^i > and pc = +00 if i^i = 0). The last equality 
in (I2.2(J|) conies from the fact that c{k,) ~ -k as k — )• 00 (see [T7] for instance), and 
the first equality comes from the correspondence (I2.17p . 

To investigate the positive solutions of equation 02.181) . we recast it into: 

which is valid as long as k 7^ 0, since c is an increasing function. This gives points (i) 
and (ii): there is no solution to 02.221) ii p < p^, and at least one solution if p > p*, 
since n 1— ?■ 4^ is a continuous function, and its infimum is p,,,. 

Let us precise now the sense of point [iii). We fix e > 0, and we suppose 
that ^ is differentiable and that for p E {pc — £,Pc) U {pc, Pc + ^); all the solutions 
of the compatibility equation satisfy {^)'{k) 7^ 0. Then, the number of solutions 
of the compatibility equation 02.221) . if finite, is odd for p e {p^ Pc + ^) and even 
for pe {pc-e,p^). 

Indeed, under these assumptions, by the intermediate value theorem, the sign 
of (-)' must be different for two successive solutions of the compatibility equa- 
tion 02.22p . Moreover, since j is unbounded (it maps its interval of definition [0, Umax) 
onto [0, +CX3)), we have 

lim 4^ = +00, (2.23) 

K-^Kmax C(k) 

SO the sign of (■^)' must be positive for the greatest solution of the compatibility 
equation 02.22p . Finally for the smallest solution, this sign must be the same as the 
sign of p — Pc- □ 

Except from these facts, since c and j are both increasing, we have no further 
direct information about this function k, i— j- j{k)/c{k). 

Remark 2.1. The results of Proposition W7^ are illustrated by FigureUi' the number 
of families of non-isotropic equilibria is given by the cardinality of the level set at p 
of the function k i— j- 4^. We see that depending on the value of p, this number 
can be zero, one, two or even more. The minimum of this function and its limiting 
value at K = provide a direct visualization of the thresholds p* and pc thanks 

to (^M>-(^M- 

We will see later on that the importance of the threshold p^ is above all due to 
a loss of stability of the uniform equilibrium, more than a change in the number of 
families of nonisotropic equilibria. And we will see that the sign of (-)'(k) which 
played a role in counting this number in the proof of point (iii) will actually play a 
stronger role to determine the stability of the nonisotropic equilibria. 

We now turn to the study of the stability of these equilibria, through the study 
of the rates of convergence. 
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Figure 1: The green, blue, red and purple curves correspond to various possible 
profiles for the function k h-s- 4^. 

3 Stability and rates of convergence to equilibria 

3.1 Main results 

We provide an overview of the most important results of this section. We emphasize 
that the results of this section are concerned with the space-homogeneous model as 
reviewed in section 12.31 and 12.41 

The first result deals with the stability of uniform equilibria. We prove that the 
critical density pc defined previously at fl2.20p acts as a threshold: 

(i) if p < Pc, then the uniform distribution is locally stable and we show that the 
solution associated to any initial distribution close enough to it converges with 
an exponential rate to the uniform distribution. 

(ii) if p > Pc then the uniform distribution is unstable, in the sense that no 
solution (except degenerate cases that we specify) can converge to the uniform 
distribution. 

The second result deals with the stability of anisotropic equilibria. As seen in 
the previous section, the anisotropic equilibria are given by the von Mises-Fisher 
distributions / = pM^^Q, defined in 02.141) . of concentration parameter k and asso- 
ciated order parameter c(k), given by the formula (12.161) . Recall that j{k) is the 
inverse function of \J\ i— ;■ A;(|J|) = T^fm- We also recall that, for a von Mises- 
Fisher distribution to be an equilibrium, the compatibility equation (12.181) i.e. the 
relation 4^ = p must be satisfied. Then: 



12 



(i) the von Mises-Fisher equilibrium is stable if (^) > where the prime denotes 
derivative with respect to k. Then, we have an exponential rate of convergence 
of the solution associated to any initial distribution close enough to one of 
the von Mises-Fisher distributions, to a (may be different) von Mises-Fisher 
distribution (with the same k, but may be different Q). 

(ii) the von Mises-Fisher equilibrium is unstable if (^j < 0. Here, the proof 
for instability relies on the fact that on any neighborhood of an unstable von 
Mises-Fisher distribution there exists a distribution which has a smaller free 
energy than the equilibrium free energy, which only depends on k, but not 
on Q. The instability follows from the time decay of the free energy. 

The main tool to prove convergence of the solution to a steady state is LaSalle's 
principle. We recall it in the next section and only sketch its proof. Indeed, the 
proof follows exactly the lines of [17] . Then, in section [331 ^6 consider stability and 
rates of convergence near uniform equilibria. Finally, in section 13. 4[ we investigate 
the same problem for non-isotropic equilibria. 

3.2 LaSalle's principle 

By the conservation relation 02.121) . we know that the free energy J-" is decreasing in 
time (and bounded from below since \J\ is bounded). LaSalle's principle states that 
the limiting value of J-" corresponds to an w-limit set of equilibria: 

Proposition 3.1. LaSalle's invariance principle: let /o be a positive measure on the 
sphere S, with mass p. We denote by J-'oo the limit of J-'{f(t)) as t -^ cxd, where f is 
the solution to the mean-field homogeneous equation (12.71) with initial condition Jq. 
Then 

(i) the set Soo = {/ G C°°(S) with mass p and s.t. V{f) = and J-'{f) = J-'oo} is 
not empty. 

(ii) f(t) converges in any H^ norm to this set of equilibria (in the following sense): 
lim dHs{f,Soo) = 0, where dH^{f,£oo) = inf \\f(t) - qWh"- 

This result has been proved in [TTj. Since the different types of equilibria are 
known, we can refine this principle to adapt it to our problem: 

Proposition 3.2. Let /o be a positive measure on the sphere S, with mass p. If no 
open interval is included in the set {n, pc{k) = j{n)}, then there exists a solution Kqo 
to the compatibility solution (I2.18P such that we have: 

lim |J/(t)| = pc{Koo) (3.1) 

and 

Vs G M, lim \\f{t) - pM,^nfit)\\H^ = 0. (3.2) 



t-^oo 
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This proposition helps us to characterize the w-hmit set by studying the sin- 
gle compatibility equation f l2.18l) . Indeed, when Koo = is the unique solution, 
Proposition 13.21 implies that / converges to the uniform distribution. Otherwise, 
two cases are possible: either Kqo = 0, and / converges to the uniform distribu- 
tion, or Koo > 0, and the w-limit set consists in the family of von Mises-Fisher 
equilibria {pM^^QjQ G §}, but the asymptotic behavior of ^f{t) is unknown. 

Proof. We first recall some useful formulas regarding functions on the sphere. Any 
function g in if* can be decomposed g = J2ege where gi is a spherical harmonic of 
degree £ (an eigenvector of — A^ for the eigenvalue £{i + n — 2), which has the form of 
a homogeneous polynomial of degree i), and this decomposition is orthogonal in H^. 
The spherical harmonics of degree 1 are the functions u ^-^ u ■ A ioT A & W"", and 
we have 

fu^udu = ^ln, i.e. yAeW, f(A-Lu)Ludu = ^A. (3.3) 

which gives that the first mode gi of g is given by the function uj \-^ nu ■ Jg, where 
the first moment Jg is defined in (12.41) . We refer to the appendix of [TTj for more 
details on these spherical harmonics. Another useful formula is 

f uVu^-A{uj)duj = - f A{uj)du, (3.4) 

where A is any tangent vector field (satisfying A{u}) ■ cj = 0). 

Since the decomposition in spherical harmonics is orthogonal in H^, we have a 
lower bound on the norm of f(t) — pMi^Q (for k ^ and fi G S) with the norm of 
its first mode: 

||/(t) - pM^nWJjs ^ [nu-iJf- J,M^,,)i-A^y[nu ■ (Jf - J.m.J] dw 

and using (13.31) . we get 

||/(t) - pM^nfns ^ n{n - iy\Jf - pci^M' (3.5) 

^ n{n - iy\\Jf\ - pc{k)\^ . (3.6) 

Since Soo consists in functions of the form pM^n with fi G § and k a solution of (I2.18p 
(and such that J-'^pM^^o) = J^oo), if we define 5*00 = {pc^n), k s.t. pc{k) = j{n)}, we 
get that the distance dH={f,£-oo) is greater than ^/n{n — iy/'^d{\Jf\, Soo), where the 
notation d{\Jf\, Joo) denotes the usual distance in M between | J/| and the set 5*00. By 
LaSalle's principle, we then have limf_^oo c?(|«/f |, •S'oo) = 0. Since |J/| is a continuous 
function, bounded in time, its limit points consist in a closed interval, which is 
included in 5*00. Obviously, if no open interval is included in the set of solutions to 
the compatibility equation (I2.18p . then no open interval is included in 6*00, and the 
limit points of | J/| are reduced to a single point pc^Hoo)- Since | J/| is bounded, this 
proves (13. ip . 

Let us now suppose that (13.20 does not hold. We can find an increasing and 
unbounded sequence t„ such that \\fitn)—pMf^^Q^.(^t„)\\H^ ^ £• By LaSalle's principle, 
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we can find gn G Soo such that ||/(t„) — gn\\ — ^ when n — )■ oo. Since gn is of the 
form pMf^^n^, we then have by the estimation (13. 6p that | J/(f„)| — pc^Kn) -^ 0, and 
so c(k„) — )• c(fi:oo)5 consequently k„ — ;■ Kqo- If «^oo 7^ 0, then we also get by fl3.5p 
that \^f{tn) - ^n\ -^ 0, so in any case, that gives that \\gn - pM^^n^.(^t„)\\H'' -^ 
(it is equal to \\pMi^,^Q^ - pM^^nj^(^t„)\\H^)- But then we obtain the convergence 
of II /(t„) - pMi^^n^(^tr,)\\H^ to 0, which is a contradiction. D 

From this proposition, the asymptotic behavior of a solution can be improved in 
two directions. First, as pointed above, the behavior of ^f{t) is unknown and we are 
left to comparing the solution to a von Mises-Fisher distribution with asymptotic 
concentration parameter n^o but local mean direction ^f{t), varying in time. If 
we are able to prove that Qf — )■ r2oo G S, then / would converge to a fixed non- 
isotropic steady-state pM^^n^o- The second improvement comes from the fact that 
Proposition 13.21 does not give information about quantitative rates of convergence 
of \Jf\ to pc(Koo), and of ||/(t) - pM^^n^(^t)\\H'' to 0, as t -> oo. 

So we now turn to the study of the behavior of the difference between the solu- 
tion / and a target equilibrium pMf^^Q^(^ty There are two tools we will use. First, a 
simple decomposition in spherical harmonics will give us an estimation in H'^ norm 
near the uniform distribution. Then we will expand the free energy J-" and its dissipa- 
tion T> around the nonisotropic target equilibrium M^^^q (^ty In case of stability, we 
will see that it gives us control on the displacement of ^f{t), allowing to get actual 
convergence to a given steady-state. We split the stability analysis into two cases: 
stability about uniform equilibrium, and stability about anisotropic equilibrium. 

3.3 Local analysis about the uniform equilibrium 

We first state the following proposition, about the instability of the uniform equi- 
librium distribution for p above the critical threshold pc- 

Proposition 3.3. Let f be a solution of (12.71) . with initial mass p. If p > Pc, and 

if Jfo 7^ 0? then we cannot have k^o = in Proposition }^.^ 

This proposition tells that the uniform equilibrium is unstable, in the sense that 
no solution of initial mass p and with a nonzero initial first moment Jj^ can converge 
to the uniform distribution. 

Proof. We first derive an estimation for the differential equation satisfied by J/ 
which will also be useful for the next proposition. 

We expand / under the form f = p + nuj-Jf + g2 {g2 consists only in spherical 
harmonics modes of degree 2 and more), and we get /§ 5^2 do; = and /§ g2 ujduj = 0. 
Let us first expand the alignment term Vtj ■ {P^±fiff) of the operator Q defined 
in (12.81) . using the fact that V^; ■ {P^±Qf) = A(^(fi/ ■ u) = —{n — l)VLf- u. We get 

V. ■ (P^xfi^/) = -p{n -l)nj-u- nV/l [(^/ ■ ^f - ^] + V. ■ (Pa.xfi/^2), (3.7) 

and we remark that the term in brackets is a spherical harmonic of degree 2, asso- 
ciated to the eigenvalue 2n of — A(^. Multiplying (12. 7p by u and integrating on the 



15 



sphere, we obtain, using (13. 71) . (13 .40 and f l3.3p (and observing that the terms /gwdw 
and Jg{u ■ QfY udu are both zero): 



d^"^^ 



n — 1 



n 



pu{\Jf\)nf + u{\Jf\)j^P^.nffduj-{n-l)r{\Jf\)Jf 
Pk{\jf\)- 



(3i 



-{n-l)r{\Jf\) 



n\J 



f\ 



Jf + ^i\Jf\) / Pu.^^fg2du. 



Using (I2.2UI) and hypothesis 12.11 we can write: 
d 



dt 



-J 



f 



(n-l)ro 1 



Pc \Jf\ 



j^Puj^92duj]jf, 



with the remainder estimation, with an appropriate constant C > 0. 

mj\)^c\j\. 



(3.9) 



(3.10) 



Equation (13.91) can be seen as ^Jj = M{t)Jf, the matrix M being a continuous 
function in time. Therefore we have uniqueness of a solution of such an equation 
(even backwards in time), and if J/^ 7^ 0, then we cannot have Jf(^t) = for t > 0. 
Now if we suppose that ||/ — p\\h' — ^ 0, then we have | J/| — ;■ and /§ Puj^g2 dw — )■ 
(as a matrix). So, for any e > 0, and for t sufficiently large, taking the dot product 
of (13. 9p with Jf, we get that 



- — I 
2dt' 



Jf\'> 



[n-l)ro{^-l]-e 

'Pc 



\J 



f\ ' 



which, for e sufficiently small, leads to an exponential growth of | J/|, and this is a 
contradiction. D 

We now turn to the study of the stability of the uniform distribution when p is 
below the critical threshold pc. We have the 

Proposition 3.4. Suppose that p < p^. We define 

A = (n-l)ro(l-^)>0. 

Pc 

Let /o be an initial condition with mass p, and f the corresponding solution of (12.71) . 
There exists S > independent of fo such that if ||/o — p\\h^ < ^, then for all t ^ 



\\f{t)-p\\Hs< 



Wfo- pWh- 
1- iWfo- pWh- 



-At 



Proof. We multiply (12.71) by {—A^)'^g2 and integrate by parts on the sphere. Us- 
ing (13.71) . (13.41) . and the fact that g2 is orthogonal to the spherical harmonics of 
degree 1, we get 



1 d 112 



-ui\Jj\)nVf\ (fi/- ^)^- 

+ iy{\Jf\) Ij^Qj ■VU-^u.y92]g2du 

-r{\Jf\) fg^i-A^r+'g^dio. 



'-A^yg2duj 
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Using the fact that the second eigenvalue of —A^ is 2n, we get 

2dt Js (3ll^ 



+ u{\Jf\) J^[Qf ■VU-^u.r92]92duj. 



We can directly compute the H'' norm of the first mode of / — p as in fl3.5p . and we 
get by orthogonal decomposition that 



11/ -P 



|2 



n(n 



-Lj K/l + 1152 



H^ 



(3.12) 



Taking the dot product of (13. 9p with n{n — 1)*J/ and summing with (13.111) . we get 
the time derivative of ||/ — p|||^s: 



Id. 

2dt 



^1^11/ - pWhs =-n{n- ir\(l - ^)\J/' 



P< 



r{\Jf\)\\92\\U. 

+ n{n - iyR{\Jf\)\Jf\^ + u{\Jf\) J^g2 % ■ Vi-A^fg^ du 

+ [n{n - 1)^ - n\2ny]u{\Jf\)\Jf\ f Qf ■ P^^^jg^doo. 

(3.13) 

Using the Poincare inequality, and again, that the second eigenvalue of — A^ is 2n, 
we get that 

IMhs+i > 2n\\g2\\ls ^{n- 1)(1 - ^)||^2||?,= . (3.14) 

Pc 

We combine the first two terms of the right-hand side of (13.131) with (I3.14p to get an 
estimation of i^||/ — p||^s in terms of a constant times ||/ — pH^'' ^^^ ^ remainder 
that we expect to be of smaller order: 



^^11/ - pWhs ^ -{n - 1)t-o(i -j)\\f- P\\l^ 



n. 



(3.15) 



where 



ns = n{n-iyR{\Jf\)\Jf\^ + u{\Jf\) f g2nfV{~A^yg2duj 
+ [n{n - ly - (2n)>(| J;|)| J^l f(Qf ■ u^g^du 



(3.16) 



[ro-r{\JMn-m-^)\\92rHs- 

Pc 



Using Lemma 2.1 of [TTj, there exists a constant Ci (independent of (72) such that 

^(72% ■ V{-A^yg2duj\ ^ CiWg^Wls. 



Together with the estimates R, v and r given by (I3.1(jp and hypothesis (12. ip . and 
the fact that the function w 1— ?■ (^2/ ■ uSy belongs to -f/^"", we can estimate every term 
of (13.160 . giving existence of constants C2, C3, such that 



7^. ^ Co 



\Jsf + \Jsn92\W + 1^/11^211^. ^ C3II/ - Pill,., 
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the last inequality coming from equation f l3.12p . Solving the differential inequal- 
ity y' ^ —\y + C^y"^ which corresponds to fl3.15p with y = ||/ — p||j?s, we get that 



y / 1/0 -\t 



A - Cay A - C-^yo 

provided that y < 5 = -^. \i yo < 5, the differential inequality ensures that y is 
decreasing and the condition y < 5 is always satisfied. In this case, we get 

which ends the proof. D 

Remark 3.1. We can indeed remove this condition of closeness of fo to p by using 
the method of JT^ in the case where p < p, where the critical threshold p is defined 
as follows: p = inf|j| ^LL (since we have c{k) ^ ^ for all k, compared to the 
definition f l2.20p -( 1^.211) of pc and p^, we see that p ^ p* ^ Pc, with a possible 
equality if for example \J\ *-^ IJI ^'^ nonincreasing) . 
We can use the special cancellation presented in JJ^: 



JVg\^ig = 0, 



where A„_i is the so-called conformal Laplacian on E>, a linear operator defined, for 
any spherical harmonic Yi of degree I, by 

K-x Ye = i{i+l)...{i + n- 2)Y,. 

Multiplying (\2.7\j by A~\(/ — p) and integrating by parts, we get the following 
conservation relation: 

Id/ n I T |2 I ||„ ||2 \_ ^/l J \\\ n /i pK\J s\) \ I T |2 , ||_ ||2 

(3.17) 

where the norms || ■ ||~_n-i and || ■ ||~_n-3 are modified Sobolev norms respectively 

equivalent to 11 ■ 11 n-i and 11 ■ 11 n-s . 
^ ''^ "h 2- '' "h~^~ 

So if p < p, equation fl3.17p can be viewed as a new entropy dissipation for the 
system, and we have global exponential convergence with rate A = (n — l)rmm(l — ^); 
where Tmin = min|j|s;pr(| J|).- 

11/ - p\\~^n^ ^ ll/o - p\\~.n^e-''\ (3.18) 

n — 1 

valid for any initial condition /q G H 2~(§) with initial mass p, whatever its dis- 
tance to p. 

Let us also remark that if p ^ p < p^,, where p* is defined in (I2.2ip . any solution 
with initial mass p converges to the uniform distribution (the unique equilibrium) , 
but we do not have an a priori global rate. We can just locally rely on Proposition \3.Ji\ 
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3.4 Local analysis about the anisotropic equilibria 

We fix K > and let p be such that k is a solution of the compatibihty equa- 
tion fl2.18p . i.e. p = 4^. In this subsection, to make notations simpler, we will not 
write the dependence on k when not necessary. 

We make an additional hypothesis on the function k: 

Hypothesis 3.1. The function \J\ i— )■ A;(| J|) is differentiable, with a derivative k' 
which is itself Lipschitz. 

We can then state a first result about the stability or instability of a non-isotropic 
solution pMf^Q, depending on the sign of (-)'. In summary, if the function k \-^ ^ 
is (non-degenerately) increasing then the corresponding equilibria are stable, while 
if it is (non-degenerately) decreasing the equilibria are unstable. For example, for 
the different cases depicted in Figure [H it is then straightforward to determine the 
stability of the different equilibria. 

Proposition 3.5. Let k > and p = 4^. We denote by T^. the value of J^^pM^o) 
(independent of Vl E Ei). 

(i) Suppose {^y{K) < 0. Then any equilibrium of the form pM^^ is unstable, 
in the following sense: in any neighborhood of pMf^Q, there exists an initial 
condition /o such that J-'{fo) < J-'k- Consequently, in that case, we cannot 
have Koo = K, in Proposition VJ/A 

(a) Suppose (-)'(k) > 0. Then the family of equilibria {pMf^^,Q G §} is stable, in 
the following sense: for all K > and s > ^^, there exists 6 > and C such 
that for all /o with mass p and with \\fo\\H^ ^ K, if ||/o — pM^n ||l2 ^ 5 for 
some f2 G S, then for all t '^ 0, we have 



\\f-pM^n,\\L-^C\\fo-pM^n 



/o 



|L2. 



Proof. We first make some preliminary computation which will also be useful for 
the following theorem. We expand the solution / of fl2.7p (with initial mass p) 
around a "moving" equilibrium pM^nj^^t)- Let us use the same notations as in [T7]: 
we write {g)M for J^g{u)M^Q^du, we denote u ■ Qf by cos 9 and we write: 

/ = M^Ufip + gi) = M^n^.{P + a{cos9 - c) + 5-2), 
where 

«= ,, ^'^'tm ■ (3.19) 

((cos6'-c)^)m 

We have {gi)M = {92) m = 0, and definition of a ensures that (w5'2)a/ = 0. The 
derivative of c with respect to k is given by 

c'{k) = (cos^ 9)m - (cos 9)Ij = ((cos 9 - cf)M. (3-20) 
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We are now ready to estimate the difference between the free energy of / and of 
the equihbrium pM^Uf ■ We have a first expansion, for the potential term of the free 
energy (12.91) : 

$(| Jjl) = <!>{pc) + k{pc)a{{cose - cf)M + k'{pc) — {{cose - cf)l + Oia^) 

J W 2 
Now, we will use the following estimation, valid for any a; G (—1, +oo): 

1(1 + x) ln(l + x) - X - ix^l ^ ||x|l (3.21) 

To get this estimation, we note that h2{x) = (1 + x) ln(l + x) — x — |x^ is such 
that /i2, /i2 3.nd /ij cancel at x = 0, and that /i2 (x) = (^^1)2 G (^l)O) for x > 0. 
Therefore Taylor's formula gives — |a;^ < h2{x) < for x > 0. For x < we have 
by the same argument h2{x) > 0, but Taylor's formula is not sufficient to have a 
uniform estimate on (—1, 0). We introduce h^ = h2 + |x'^. By induction from 2 = 3 

to 2 = 1 we have that h]^ as a unique root 7j in (—1,0), with 73 > 72 > 7i. 
Since h'^{x) — )■ — 00 as x — )■ —1, h^ is decreasing on (—1, 71) and increasing on (71, 0), 
but we have h-si^—l) = /i3(0) = so /13 < on (—1,0), which ends the derivation 

Using (13.211) with x = — , we have that 



J§ 



flnfdu ={{p + gi)[\n{l + ^) + \n{pM^n,)])M 

p 

= {pHpM^nf))M + {ticosegi)M + y(91)m + Oi{\9if)M) 

pM^Uf \n{pM^n^.)duj + aKc' + ^[«^c' + (^2^)^/] + 0{{\gi\'^)M). 
Finally we get 

= ^[$(;)'«' + i9l)M] + 0{{\g,\')M). (3.22) 

Now, we prove (i). We simply take a sufficiently small and g2 = 0, and the estima- 
tion (13.221) gives the result. Indeed, since c and j are increasing functions of k, the 
leading order coefficient in (13.221) . which is Y^ii) ' ^^ negative by the assumption. 
We now turn to point (ii). We will use the following simple lemma, the proof of 
which is left to the reader. 

Lemma 1. Suppose x{t) ^ is a continuous function and y{t) is a decreasing 
function satisfying 

\x{t)-y{t)\ ^Cx(t)i+",Vt^O, 
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for some positive constants C and e. Then there exist 5 > and C such that, 
if x{0) ^ S, then 

y{t) ^ 0, and \x{t) - y{t)\ ^ Cy{tY^' ,\/t ^ 0. 

By Sobolev embedding, Sobolev interpolation, and the uniform bounds of The- 
orem [H we have 

ll^illoo ^ C\\gi\\ n_^ ^ CWg^WH^hifL^ ^ Ci{{gl)Mr. (3.23) 

for some e > 0, and where the constant C\ depends only on K (the constant in 
the statement of the proposition, which is an upper bound for ||/o||_ff05 ^i ^ ^'^^ 
the coefficients v and r of the model. We will denote by Ci such a constant in the 
following of the proof. 

We define x{t) = j^iyH)' a"^ + (^Da/] and y{t) = T{f)-T^. Together with the 
estimate (I3.22p . since {gl)M = c'a^ + (s'Dm, and {^Y > 0, we can apply Lemma [H 
It gives us that if {g'Du is initially sufficiently small, then J^{f) ^ T^ and we have 

<t) = ^[^(;)' «^ + {qDm] = Hf) - -^. + o{{Hf) - j^^y^n- 

z,p J c 

Now, using the fact that x(t), {g'Du and ||/ — pM^^j Ili2 are equivalent quantities 
(up to a multiplicative constant) and the estimate (I3.22J) . we get that 

\\f - pM^^^Wl, ^C^xit) ^a.iJ'U) - ^>.)- (3.24) 

Using the fact that J^{f) — T^, is decreasing in time, and the same equivalent quan- 
tities, we finally get 

\\f-pM^^^\l, ^ C'3(^(/o) - J-.) < C4II/0 -pM,n,J|i2. 

This completes the proof, with the simple remark that, as in the proof of proposi- 
tion [2121 we can control 1^2 — fi/gl by \fQ — pM^£i\i2 (using the formula (13.51) ). Then 
we can also control the quantities \p{MnQ. — ^kQ.j )IU2 and ||/o — pM^Q^. \\l2, and 
finally the initial value of {g'Djvh by this quantity ||/o — pM^^W^-z. D 

We can now turn to the study of the rate of convergence to equilibria when it 
is stable (in the case (-)' > 0). The main result is the following theorem, which 
also gives a stronger stability result, in any Sobolev space H^ with s > ^^. Let 
us remark that this theorem is an improvement compared to the results of [T7], in 
the case where r is constant and z^(| J|) is proportional to \J\. In what follows, we 
call constant a quantity which does not depend on the initial condition /o (that is 
to say, it depends only on s, k, n and the coefficients of the equation u and r). 

Theorem 2. Suppose (^)'(/«) > 0. Then, for all s > ^^, there exist constants S > 
and C > such that for any /o with mass p satisfying ||/o — pM^nUn^ < S for 
some Q eS, there exists Q^o G § such that 

11/ - pM^njHs ^ C\\fo - pM^nllH^e-^', 
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where the rate is given by 

A = ^A.(^)'. (3.25) 

J c 

The constant A^ is the best constant for the following weighted Poincare inequality 
(see the appendix of fTTy for more details on this constant, which does not depend 
on VL): 

{\^..g?)M>K{{g-{g)Mf)M. (3.26) 

We first outline the key steps. Firstly, we want to get a lower bound for the dissipa- 
tion term 'D{f) in terms of T{f) — J-'k, in order to get a Gronwall inequality coming 
from the conservation relation (12.121) . After a few computations, we get 

!?(/) ^ 2A(^(/) - J-.) + 0((^(/) - :F.f^^). 

With this lower bound, we obtain exponential decay of T{f) — J-'k (with rate 2A), 
which also gives exponential decay of ||/ — M^j^^, 11^2 (with rate A) in virtue of ([32 



We also prove that we can control the displacement (if by J{g2)M- Hence we get 
that Qf is also converging exponentially fast towards some Qqo ^ S (with the same 
rate A). After linearizing the kinetic equation (12.71) around this equilibrium pMi^n^, 
an energy estimate for a norm equivalent to the H^ norm gives then the exponential 
convergence for ||/ — Msf^^Hz/s with the same rate A. 

We now give the detailed proof. 

Proof of Theorem\^ We fix s > ^^ and we suppose (^)'(k) > 0. We recall the 
notations of the proof of Proposition 13. 5t 

/ = M«n^(p + 5fi) = M«,n^(p + a (cos 6^ - c) +5-2), 

where cos^ = 00 -Vlf and a, defined in (I3.19p . is such that 

I J/l = pc + a{{cos9 — cY)m = j + cue , (3.27) 

thanks to (l3.2Up . We have that {gi)M = (fi'2)A/ = 0, and (^(72)1/ = 0. 
The proof will be divided in three propositions. 



Proposition 3.6. There exist constants S > 0, e > and C such that, if initially, 
we have {g\)M < ^ (md ||/o — M^^q Wh" ^ 1, then for all time, we have 

V{f) ^ 2A(^(/) - J-.) - C(^(/) - J•.)l+^ 

where the rate is given by (I3.25p .' A = ^-^r^Ai^{^)' . 

Proof. We apply the stability results of the second part of Proposition 13.51 with the 
constant K being 1 + UpAf^nj Wh" (this does not depend on fif^). This gives us 



22 



constants Si > 0, e > 0, Ci, C2 such that if we have initially {gf)M < ^i, then (see 
formulas (E23D-(E21) 



\gi\\oo^Ci{glY^, (3.28) 



1 rCc'.ix/ 



o}^0'^ +^9i)M] - {Hf)-^.) ^ C.iTin-J'.r^, (3.29) 

Zp J c 

{gDM^Csinn-J",). (3.30) 

We get, using the definition 02. lip : 

Vif) = Ti\Jf\){ip + gi)\V4Hp + gi)-iki\Jj\)-K)u-nj]\')M 

= ri\Jf\){^—\V^gi\' + (p + ^i)(A:(l^/l) - «:)'| V.(c^ ■ nf)\')M 

p + gi 

- 2ri\Jf\){V^gi ■ iki\Jj\) - k)V^{u ■ %))m. 
Using the fact that — p— ^ ^{p — \\gi\\oo), ^^ obtain 



1 



^(/) > ri\jf\)ip - \\g,\\^){ \v^gi\' + (Ml J/l) - f^rivu^ ■ nj)\') 



M 



p^ 

- 2r(| J;|)(V^(7i • (A;(| J;|) - ^)V^{uj ■ Qf))M 

VU) > ^(l>//l)(P - ll^l||oo)(|V.[| - (Ml J/l) - ^)^ ■ ^^M 

+ ^(l^/l)-|l^illoo(A;(| J/l) - K){gi{K\VUo^ ■ fi/)|2 - (n - l)iu ■ fi/))M. 

where we used Green's formula to evaluate (V^(7i ■ V^(a; ■ Qf))M- 

First of all, using the definition (I3.27P we can get that |fc(| J/|) — k — a%| ^ 6*40; 
for a constant C4. Then we use the Poincare inequality 03.260 : 

(|V^^P)a/^A«(((7-(^)m)')m. 



2 



) 



Hence, since |q;| is controlled by JigDhi (we recall that {gf)M = c'a'^ + ((^Dm), and 
since we also have ||J/| — j\ ^ Csl^l for a constant C5, we get 

V{f) ^ AM\Jf\)iP - \\9i\U{\- - {k{\Jf\) - «:)(cose - c)\')m 



-C'6||5'i||oo(5'?) 



M 



^ A.r(j)p(|^ + «(- - ^)(cos0 - c)\')m - Cr\\g4oo{gl)M 
P P 3 

A-kT(JJ |. C C ^^J^^/^^2 2 , / 2\ 1 r~< U II / 2\ 

>M, 



P 



[jjjyX-yy^ +(9i)M]-Cr\\g,\Ugf). 



where Cq and Cj are constants. Together with the fact that 4('^)' ^ 1 (this is 
equivalent to jc' ^ 0), and with equations O3.28j) - O3.30p . this ends the proof. D 
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Proposition 3.7. There exist positive constants C,C and 6 such that if initially, 
we have {qDm < ^ and ||/o — pM^Q^. \\h<> ^ 1, then for all time, we have 

11/ - pM^n, \\l^ ^ C\\fo - pM^n,^^ \\l2 e"^*, 

and furthermore, there exists Q^o ^ S such that for all time, we have 

\Qf-Q^\^C\\fo-pM^njL^e-^K (3.31) 

Proof. By Proposition 13.6^ using the expression {gf)M = c'a^ + {02) m and inequali- 
ties (13.291) and (13.301) . we get that there exist constants Si > and Ci, C2, and 6*2 > 
such that if {gf)M < ^1, then J-'{f) ^ J-'k, and for all time, 

^(^(/) - J-.) = -V{f) ^ -2A(^(/) - J-.) + Ci(^(/) - J•«)l+^ (3.32) 

C2iJ'if) - J".) ^ {qDm ^ C2{Hf) - ^n). (3.33) 

Solving the differential inequality (13.321) for F{f)—J^K sufficiently small, we get that, 
up to taking 62 < 61, if {qDm < ^2, we get a constant C3 such that 



^(/)-J-«^C3(^(/o)-J'.)e 



^2 At 



This gives the first part of the proposition, with (I3.33p . and the fact that there exists 
constants C4, and C4 such that 



C4II/ -pA4n,||L^ ^ ^{qDm ^ C4f - pM^njh^. 
Now we compute the time derivative of Qf, using ^^f = rrf-Pfi-L^'^/ and (13. 8|) : 

4^/ = ^^^Pn|(Ofi/(p + «(cos^-c) +^2))m = -^^^Pnf{cos9ug2). 
So there exist constants C5 and Cq such that 



which, after integration in time, gives the second part of the proposition. D 

We can now prove the last step which leads to Theorem [21 

Proposition 3.8. There exist constants 5 > and C > 0, such that for any initial 
condition /o with mass p satisfying ||/o — pM^Q^. \\h<' < S, there exists r2oo G § such 
that 

\\f-pM,nJ\Hs^C\\fo-pM^njH^e-^'. 

Proof. All along this proof we will use the symbol x to denote quantities of the 
same order: for a and b two nonnegative quantities, a ^ h means that there exist 
two positive constants Ci, C2 such that Cia ^ 6 ^ 6*20. 

By the estimation {g\)M ^ 11/ — P^kQ/ ||i2 (since the weight M^q is bounded 
above and below), and by a simple Sobolev embedding (L^ C H^), there exists a 
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constant Si > such that if ||/o — pMi^n^, \\ho < ^i, then we are in the hypotheses 
of Proposition 13. 71 We suppose we are in that case and we can then go back to 
the original equation and perform a hnear analysis around pMi^Q^. We will now 
write {g)M for J^gM^n^du. 

If we write f = {p + g)Mf^Q^, then the equation becomes 

dtg = -r{\Jf\)Lg - A{t) ■ V^g + B{oj) ■ A{t){p + g), (3.34) 

where 

Lg = -— V^ ■ (M^n^V^g) = -{A^g + nVt^ ■ V^g), 

Ait) = u{\Js\)^f - T{\Jf\)K^^, 

B{uj) = [n — l)u} — kP^±Qoo- 

Let us remark that the linear operator L is a coercive self adjoint operator for the 
inner product {gi,g2) H- {gig2)M (also denoted {gi,g2)M in the following), on the 
space LIj C L^ of functions g such that {g)M = (thanks to the Poincare inequal- 
ity (I3.26P ). Indeed we have 

{gi,Lg2)M = (V^s-i ■ V^g2) hi- 
lt is classical to prove that the inverse of L is a positive selfadjoint compact operator 
of L\[. Hence, by spectral decomposition, we can define the operator L^ , and use it 
to define a new Sobolev norm by 



\9\?Hs = {g,L'g)M. 



We will use a lemma (the proof of which is postponed at the end of this section) 
about estimations for this norm, and about a commutator estimate: 

Lemma 2. For s ^ 0, we have \\g\\H= ^ lls'll^"? for functions g in Hlj = H^ n L\j. 
Furthermore, forg G Hl^, the (vector valued) quantity {U gV u)g) m is well defined 
and there is a constant C such we have: 

\{L'gV^g)M\^C\\g\\%. (3.35) 

M 

We will also use the following Poincare estimate, for g G -ff|^, with the same 
constant A^ as in (I3.26P : 



{g.L^+'g)M = (|V(Lf^)p)M ^ K{{L'^gf)M = ^M\% 



We now multiply the equation 03.341) by L'^g and integrate with respect to the 
measure M,^Q_^duj. We get 



1 d 



^Ilk^-^ll-^/DA.II^IlL +1^(^)1(^111(711^. +||(7||H^,l|5H(P + ^)bjJ, 



where \\B{uj){p + g)\\Hs denotes the maximum of ||e ■ B{uj){p + g)\\Hs for e G S. 
Since w i— ;■ e-B{u) is smooth, the multiplication by e-B{(jj) is a continuous operator 
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from HIj to H'^ when s is an integer, so by interpolation this is true for all s. 
Therefore, we get a constant C2 such that for all g G Hlj, we have 

\\B{uj)gUs^^^C2\\g\\i,.^^. (3.36) 

We finally get 

JtMtii, ^ -r{\Jf\)K\\g\\Hi^ + I^WI((C'i + C2) hWiii^ + ||i?(a;)p||^.;. 
Now, applying Proposition 13.71 there exist constants C^.C^.C^ such that 
\A{t)\ ^ v{\Jfmj - fiool + H\Jf\) - r{\Jf\)K]\n^\ 

^ U{\jf\)\^f - ^^\ + T{\Jf\)[k{\Jf\) - kU{K))] 

^ Cal^ - f^ool + C4II J/l - j(«^)| 

in virtue of (13.61) . By the same argument, we get, for a constant Ce, that 

-\t 



\r{\J}\)-r{j)\^C4U-pM^nA\ 



e 



so we finally obtain, together with a uniform bound on 1 1(711^/8 coming from The- 
orem [T] (and independent of /o since ||/o — pM^Q^ \\h^ < (^i), a constant C7 such 
that 

We solve this inequality and we get 

ll^lli/f, ^ II^oIIh; exp(-r(j)A^t) + C7||/o-pM«n,J|i/?,^ ^, 



'm r(j)A.-A 



and this gives the final estimation, using the fact that A < r(j)AK (this is equivalent, 
by definition 03.251) of A, to (-)' < — , and we indeed have j c' > 0), and that 



^\\fo-pALnjH^^^+Cs\nj,-a 



'-f^M 



^Cgll/o-pAf, 



KQf^WH^ 



in virtue of (13.61) . Lemma El and (I3.3ip (we have s > ^^ so L^ C H'^ is a continuous 
embedding) . D 

Finally, Proposition 13.81 can be refined since, thanks to the estimation (13.61) . 
we only need to control ||/o — pMi^^W^s for a given fi G S in order to ensure 
that Wfo — pMi^Q ll^s is sufficiently small, and this ends the proof of Theorem [2J D 

Proof of Lemma [3 We first define the space i/fj as the completion of C°° (S) fl L\.j 
fo^ II ■ II ij" • The first estimate (which amounts to prove that Hfi'll^fs x ||5'||//s 
for smooth functions g G C°°(§) fl L\j) is true when s is an integer: indeed L* 
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and (— A^)'' are simple differential operators (of degree 2s), and these estimates can 
be done by induction on s: when s = 2p is even, we write 



'{g,L^g)M = \\L^g\\l.><\\L^g\\l 



\9\\h 



-AF^I 



-Afg 



|2 



In the first case, L is decomposed as (— A^^) — kQod ■ V^; to estimate ||L^(7||2 in 
terms of 11(711/^8, and in the second case — A^; is decomposed as L + k^oo ■ V,j to 
estimate ||(— A)^5f||^2 in terms of {g^L^g). When s = 2p + 1 is odd, the same 
argument applies, writing 



M 



Mis 



IV. 



|V.(L^^)| 



f 2 



|V.(LP^)| 






Finally, the general case is done by interpolation, for s = n + 6, with 6 G (0, 1). 
We refer the reader to [2S] for an introduction to interpolation spaces, and we will 
denote (Fi, F2)(5ip) the interpolation space between Fi and F2 using the real inter- 
polation method. Using the so-called /f-method (see [251 Lecture 22]), it consists 
in the space of elements u E Fi + F2 such that ||m||0,p < +00, together with the 



norm 



Mid,; 



iP' 



where 

t / 



with K(t, u) 



inf 

UlGFi,U2eF2 



I^iIIfi +^||m2||f2- 



We will use the following result (see [2H1 Lemma 23.1]): if (X, fj.) is a measured space 
and Wq, Wi are two weight functions, we have 



{L\wodfi),L\widfi))^g^2) = L\wyw',dfi), 



(3.37) 



where, for a weight function w ^ 0, the weighted space L^(wdyu) denotes the func- 

i2(t«dM) = Jxh'^{x)w{x)d^{x) is finite. Now if {gi)i(zN is an 



tions h such that 

orthonormal basis (for the dot product (■, ■)m) of eigenvectors of L (associated to 
the eigenvalues (Aj)), it is easy to see that the map h \-^ {{h,gi)M)ieN is an isom- 
etry between H^j and the weighted £^ space with weight (Aj)jgN (it corresponds 
to L'^{wdji) where X = N, /i is the counting measure, and w{i) = Aj). There- 
fore, we obtain with (I3.37P that Hlj = {H^, H2^^)(e,2), and by the same argu- 



ment H' = (i7", H 



n+U 



)(e,2) 



. So we finally get, with equivalence of norms: 



H' n Ll, 



(//" n lIj, /7"+i n L^ 



m){B,2) 



(^M'^A/^)(e,2) = H 



M- 



To get the estimation (13.351) . we first observe that it is a commutator estimate. 
Indeed, by integration by parts for a given e G S, we get that the adjoint operator 
of e- Vo; (for (-, ■)m) is — e- Va; + e-i?(a;), where B{u) = {n — l)u — KPi^±flao (the same 
expression as in the proof of Proposition 13. 8p . So, splitting the left part of (I3.35P in 
two halves, we are led to show that for g G Hl^, we have 



-\{g[L^,V^]g)M 



{B{uj)gV 



)M 



^C{g,V 



Hi- 
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Using fl3.36p . it is equivalent to find a constant C such that for all g G -ff|^, we have 

\{g[L\Vj^g)M\^C{g,Vg)M. (3.38) 

In the case s = 1, by using Schwartz Theorem, we see that [L, V^;] = [— A^, V^^]. It 
is proven in Lemma 2.1 of |T7j that f l3.35p is true in the limit case where k = 0. This 
means that [(— A^)"*, V^;] is an operator of degree 2s. In particular [— A^, Vcj] is a 
differential operator of degree 2. Actually, using Lemma A.b of [TTj, it is possible to 
get that 

[-A^,Vj=2c^A^-(n-3)V^. 

This directly gives the estimate (13.381) when s = 1. We obtain the estimate when s 
is an integer with the formula [Lp+\ V^;] = Eg=o^^"^[^, V^lL''. 

The proof in the general case relies on a resolvent formula for the operator A^ , 
when 9 belongs to (0, 1), and A : 'D{A) G H -^ H is a strictly positive operator of 
a Hilbert space H with a complete basis of eigenvectors (see [221 Remark V-3.50]): 



A^ = / t^(r^ - (t + A)-i)dt. 

TT JO 



This formula can be checked on an orthonormal basis of eigenvectors of A, and relies 
on the fact that, for A > 0, we have 



'') 



oo 



f -]dt = X^ 



U t + X^ Jo 1 + t 

The fact that this last integral is equal to . "" „ for < 6* < 1 is classical, and can 

<= ^ sm Try ' 

be done by the method of residues. 

We then have, for another operator B (with dense domain for {t + A)^^B 
and B{t + A)'^ for t > 0) 

[A\B] = '-^^ rt%B,it + Ar']dt 

TX Jo 

= ^^^ rt\t+Ar'[A,B]{t+Ar'dt. 

TT Jo 

We can apply this result to A = L™ with H = L\j, B = e ■ V(^ for a fixed e G S, 
and ^ = — forO<s<m and we get, using the fact that {t + L™)~^ is self- 
adjoint in H (and bounded, so all smooth functions are in the domain of {t + A)^^ B 
and B{t + A)-^), 



siUTT^ 

sin7r6' 



sin TTr? r^^ 

\{g[L^,e-W^]g)M\^ / t^K^ (t + L™)-i[L-, e ■ V.](t + L™)-^!/^^ 

vr JO 

/•oo 

Jo 



^ _^ / f\{(t + L"'r'g[L"',e-V^]{t + L"'r'g)M\dt 

71 



sin it9 r°° 

n 
sin 7i9 roo 



n Jo 





)0 

m\— Irm/j. , rm\— 1, 



M 



dt 



^ Cm / f{g, {t + L'")-^L™(t + L"')-'g)Mdt. 

IT Jo 
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But as before, it 


is easy to see that 




io \t + Xy Jo {l + ty~ smnO 


and then 


9A' = '''' ""^ r t\t + A)-'A{t + A)-Mt. 

71 Jo 


Finally, we get 






m 



wich ends the proof of Lemma [2] □ 

4 Phase transitions 



4.1 Application of the previous theory to two special cases 

In the previous section, we have stated results which are valid for all possible be 

c(k) 



haviors of the function k i— j- 4^. In particular, the number of branches of equilibria 



can be arbitrary. 

In this section, we apply the previous theory to two typical examples: 

(i) The function k i— )• 4^ is increasing. In this case, there exists only one branch of 
stable von Mises-Fisher equilibria. The uniform equilibria are stable for p < pc, 
where pc = limK-5.0 4fl, and become unstable for p > pc- The von Mises-Fisher 
equilibria only exist for p > Pc and are stable. This corresponds to a second- 
order phase transition. We will provide details and a determination of the 
critical exponent of this phase transition in section 14. 2[ 

(ii) The function k h- 4^ is unimodal, i.e. there exists k^ such that this function is 
decreasing on [0, n^] and increasing on [k^, 00). Then, another critical density 
is defined by p* = 4^- Then we have the following situation: 

a. if p G (p^,, Pc), there exist two branches of von Mises-Fisher equilibria, and 
therefore, three types of equilibria if we include the uniform distribution. 
Both the uniform distribution and the von Mises-Fisher distribution with 
the largest k are stable while the von Mises-Fisher distribution with 
intermediate k is unstable. 

b. if p < p*, there exists only one equilibrium, the uniform one, which is 
stable. 

c. if p > Pc there exist two types of equilibria, the uniform equilibrium 
which is unstable and the von Mises-Fisher equilibria which are stable. 

This situation corresponds to a first-order phase transition and is depicted in 
section 14.31 where phase diagrams for both the two-dimensional and three- 
dimensional cases are given. The major feature of first-order phase transitions 
is the phenomenon of hysteresis, which will be illustrated by numerical simu- 
lations in section 14.31 

For references to phase transitions, we refer the reader to 
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4.2 Second order phase transition 

Let us now focus on the case where we always have (-)' > for all k > (this 
corresponds for example to the upper two curves of Figure [T]). In this case, the 
compatibility equation fl2.22p has a unique positive solution for p > pc- With the 
results of the previous subsection about stability and rates of convergence, we obtain 
the behavior of the solution. 

Proposition 4.1. Let /o be an initial condition with mass p, and f the corresponding 
solution of (12.71) . We suppose that (^)' > for all k > 0. Then: 

(i) If p < Pc then f converges exponentially fast towards the uniform distribu- 
tion foo = P- 

(a) If p = Pc, then f converges to the uniform distribution foo = P- 

(Hi) If p > Pc and Jf^ ^ 0, then there exists Vtao such that f converges exponentially 
fast to the von Mises-Fisher distribution f^o = pMi^^^, where k, > is the 
unique positive solution to the equation pc^n) = j{k). 

Proof. This is a direct application of Propositions 13.2113.41 and Theorem [21 □ 

Remark 4.1. 

(i) When p > pc, the special case where Jfg = leads to the study of heat equa- 
tion dtf = tqA^/. Its solution converges exponentially fast to the uniform 
distribution, but this solution is not stable under small perturbation of the ini- 
tial condition. 

(a) For some particular choice of the coefficients, as in fTTf - it is also possible to 
get a polynomial rate of convergence in the second case p = pc. For example 
when j{k) = k, we have \\f — p\\ ^ Ct^2 for t sufficiently large. 

We now describe the phase transition phenomena by studying the order param- 
eter of the asymptotic equilibrium a = ■^°° , as a function of the initial density p. 

We have a(p) = if p ^ pc, and a(p) = c(k) for p > pc, where k > is the 
unique positive solution to the equation pc{K,) = j{n). This is a positive continuous 
increasing function for p > pc- This is usually described as a continuous phase 
transition, also called second order phase transition. 

Definition 4.1. We say that /3 is the critical exponent of the phase transition if 
there exists Oq > such that 

a{p) ~ ao{p - pcY, as pA- pc- 

This critical exponent /3 can take arbitrary values in (0,1], as can be seen by 
taking k such that j{k) = c{k,){1 + k^) (we recall that k is the inverse function of j, 
see Hypothesis 12. 2p . Indeed in this case, the function k is well defined (its inverse j is 
increasing), and satisfies Hypothesis 13.11 (if B ^ 1). We then have (^)' = 4k^~ > 0, 
and the conclusions of Proposition 14. II apply, with pc = 1. Finally, the compatibility 
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equation pc{K.) = i{n) becomes p= (I + k/^), i.e. k = {p — 1)^. And since c[k) ~ -k 
when K — )■ 0, we get: 



a(p) = ciip - lY) ~ -ip - 1)'^ as p 4 1 
^ ^ n 



More generally, we can give the expression of the critical exponent in terms of the 
expansion of k in the neighborhood of 0. 

Proposition 4.2. We suppose, as in Proposition \4.1\ that (-)' > for all n > 0. 

We assume an expansion of k is given under the following form: 

^Ml = !L-a\J\'^ + o(| J|^) as\J\^ 0, (4.1) 

\J\ Pc 

with q ^ 1 (see Hypothesis \3.1\) and a ^W. 

(i) If q < 2 and a j^ 0, then a > and we have a critical exponent given by (3 = -. 
(a) If q > 2, the critical exponent is given by /3 = ^. 

2 2 

(Hi) If q = 2 and a ^ — 3 ." ^. , then a > — 3 " , and the critical exponent is given 
by P = \- In the special case where 

^ = - + V^i^i' - ^^1^1' + ^d^i') '' 1^1 ^ 0' 

\J\ pc Pt{n + ^) 

with 2 < p < 4 and 02 7^ 0, then 02 > and we have a critical exponent given 

It is also possible to give more precise conditions for a higher order expansion 
of k in order to have a critical exponent less than |, the point (iii) of this proposition 
is just an example of how to get an exponent less than |. We will only detail the 
proofs of the first two points, the last one can be done in the same way, with more 
computations, which are left to the reader. 

Proof. We only detail the first two points, the last one is done in the same way, with 
more complicate computations. We recall that k{j{K)) = k by definition of j. So we 
get that K ~ "^lilnl ag ;^ _). Q^ And using fl4.ll) . we obtain 

_K_ ^ k{j{K)) 

Furthermore, we have ^^ — n~ n'^(n+2) ^^ ~^ 0{k^) (see [IT], Remark 3.5) as k — )■ 0. 
So we get, as K — )■ 0: 

p ]{k) j{k) k Pc\ ^''' "(^^+2) ; V ) 
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n 


fl^Pc 

- a 




Pc 


V n , 


/ ^ 



So since k i— ;■ 4^ is decreasing, if g < 2 and a 7^ we have a > 0. In this case, 
we get that p = pc{l + a{^)'^~^^ '^^) + o(/t'^) as k — ;• 0. Hence, as p — t- pc, we 
have K, ~ -^ — TTrip ~ Pc)''- Since c(k) ~ - as k — )• 0, we obtain (i). 

For the same reason, if g > 2, we get p = pc{l + n(n+2) ^'^) + o(k™™^'^'^)) as k — ;• 0, 



n{n+2) 



and then k ~ J^^ — '-(p — p^) as k — )• 0, which proves point (ii). D 

The hypothesis in Proposition 14.11 is not exphcit in terms of the ahgnment and 
diffusion rates u and r. We have a more direct criterion in terms of k which is given 
below (but which is more restricted in terms of the critical exponents that can be 
attained) . 



Lemma 3. // , ,, is a non-increasing function of \J\, then we have (^)' > /( 



1^1 



or 



all n > 0. In this case, the critical exponent, if it exists, can only take values in [|, 1]. 
Proof We have that £(^) < for k > (see ^7\). Then 



ck{j)J c\k{j)J \cj k{j) 

since {-J^Y ^ (j is an increasing function of k and -|j| is a non-increasing 
function of | J|). Now if we suppose that there is a critical exponent /3 according 
to Definition 14.11 we get, using the fact that a{p) = c{k) ~ - as k — )■ 0, that 

- = - — an^ + o(k^) as k — )• 0, with a = {pc)~'^{nao)~^ . We then have 



^ = - - = (n + -^k' + 0(«:^))(- - a J + o0)) 

Pc (n + 2)pe 
Then B cannot be less than ^, otherwise the function ^■'iJ' could not be a nonin- 

f^ 2' j(k) 

creasing function of k in the neighborhood of 0. D 

Remark 4.2. When this criterion is satisfied, we can also use the result of Re- 
mark \y.l[ Indeed, in that case we get easily that p = pc, and we obtain that 
for any p < pc, there is a global rate of decay for the modified H^^2- norm: for 
all /o G H~^~{'B), we have the estimation 03.18p . 



4.3 First order phase transition and hysteresis 

We now turn to a specific example, where all the features presented in the stability 
study can be seen. We focus on the case where i/(| J|) = | J|, as in [T7|, but we now 
take r(| J|) = 1/(1 + \J\). From the modeling point of view, this can be related to 
the Vicsek model with vectorial noise (also called extrinsic noise) [Il[3], since in that 
case the intensity of the effective noise is decreasing when the neighbors are well 
aligned. 
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In this case, we have /i;(|</|) = \J\ + | Jp, so the assumptions of Lemma E] are not 
fulfilled, and the function j is given by j{k) = |(-\/l + 4:K — 1). 
Expanding ^ when k is large or k is close to 0, we get 



c \^ + 0{l) 



as K — )• 0, 

as K — > oo. 



Consequently, there exist more than one family of non-isotropic equilibria when p is 
close to pc = n (and p < pc). 

The function k h- 4^ can be computed numerically. The results are displayed 
in Figure E] in dimensions n = 2 and n = 3. 

4 




Figure 2: The function k i— ?■ 4^, in dimensions 2 and 3. 

We observe the following features: 

• There exists a unique critical point n^, for the function -, corresponding to 
its global minimum p^, (we obtain numerically p* ~ 1.3726 and n^, ~ 1.2619 
if n = 2, and p* ^ 1.8602 and k, ^ 1.9014 if ra = 3). 

• The function - is strictly decreasing in [0, «=„) and strictly increasing in (k*, oo). 

We conjecture that this is the exact behavior of the function -, called unimodality. 
From these properties, it follows that the solution associated to an initial condition /o 
with mass p can exhibit different types of behavior, depending on the three following 
regimes for p. 

Proposition 4.3. We assume that the function - is unimodal, as described above. 
Then we have the following hysteresis phenomenon: 
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(i) If p < p*, then the solution converges exponentially fast to the uniform equi- 
librium foo = P- 

(a) If p* < p < n, there are two families of stable solutions: either the uniform 
equilibrium f = p or the von Mises-Fisher distributions of the form pM^^ 
where k, is the unique solution with n > n^, of the compatibility equation f l2.18p 
and fi G S. // /o is sufficiently close to one of these equilibria, there is expo- 
nential convergence to an equilibrium of the same family. 

The von Mises-Fisher distributions of the other family (corresponding to so- 
lution of (12.181) such that < k < k^,) are unstable in the sense given in 
Proposition \3.5\ 

(Hi) If p > n and Jf^ ^ 0, then there exists Q^o ^ S such that f converges exponen- 
tially fast to the von Mises-Fisher distribution pM,^^^, where n is the unique 
positive solution to the compatibility equation pc{k) = j{k). 

Proof. Again, it is a direct application of Propositions I3.2tl3.4l and Theorem El □ 

Remark 4.3. 

(i) At the critical point p = p^, the uniform equilibrium is stable (and for any 
initial condition sufficiently close to it, the solution converges exponentially fast 
to it), but the stability of the family of von Mises-Fisher distributions p^^Mi^^Q, 
for fi G S, is unknown. 

(a) At the critical point p = n, the family of von Mises-Fisher distributions of the 
form nMf^^Q, for f2 G S and where n^ is the unique positive solution of fl2.18|] . 
is stable. For any initial condition sufficiently close to uM^^q for some fi G §, 
there exists f2oo such that the solution converges exponentially fast to nMi^^n^ . 
However, in this case, the stability of the uniform distribution f = n is un- 
known. 

(Hi) As previously, in the special case Jf^ = 0, the equation reduces to the heat 
equation and the solution converges to the uniform equilibrium. 

The order parameter ci as a function of p (i.e. Ci(p) = c{k) with p = |^) is 
depicted in Figure E] for dimension 2 or 3. The dashed hues corresponds to branches 
of equihbria which are unstable. 

The hysteresis phenomenon can be described by the hysteresis loop. If the pa- 
rameter p starts from a value less than p^,, and increases slowly, the only stable distri- 
bution is initially the uniform distribution and it remains stable, until p reaches the 
critical value pc- For p > pc, the only stable equilibria are the von Mises-Fisher dis- 
tributions. The order parameter then jumps from to Ci(pc)- If then the density p 
is further decreased slowly, the von Mises-Fisher distributions are stable until p 
reaches p^, back. For p < p*, the only stable equilibrium is the uniform distribution, 
and the order parameter jumps from Ci(p=k) to 0. The order parameter spans an 
oriented loop called hysteresis loop. 

This hysteresis loop can be observed numerically at the kinetic level or at the 
particle level. The plots of the order parameter for such numerical simulations are 
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Figure 3: Phase diagram of the model with hysteresis, in dimensions 2 and 3. 

given by Figures H] and O The details of the numerical simulations are provided 
in appendix [Al The key point to be able to perform these numerical simulations 
while varying the parameter p in time is to rescale the equation in order to see the 
parameter p as a coefficient of this new equation, and not anymore as the mass of 
the initial condition (normalized to be a probability measure). 

We can also obtain the theoretical diagrams for the free energy and the rates 
of convergences. For this particular example, the free energies J^{p) and J-'^ (cor- 
responding respectively to the uniform distribution and to the von Mises-Fisher 
distribution pM^Q for a positive solution k of the compatibility equation (I2.18P ) are 
given by 



J^{p) = plnp, 
j; = plnp + {p\nM^n)M - -f 

= plnp -pin /e"^°^Mu;-- 
J 6 



-3^ 



7J«^- 



The plots of these functions are depicted in dimensions 2 and 3 in the left plot 
of Figure El Since the functions are very close in some range, we magnify the 
difference J-'^ — J-'{p) in the right plot of Figure El The dashed lines correspond to 
unstable branches of equilibria. 

We observe that the free energy of the unstable non-isotropic equilibria (in 
dashed line) is always greater than the one of the uniform distribution. There ex- 
ist pi e {p*,pc) and a corresponding solution ki of the compatibility solution 02.180 
(with Ki > /€*, corresponding to a stable family of non-isotropic equilibria) such 
that J-'k^ = J-'(pi). If p < pi, the global minimizer of the free energy is the uniform 
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Figure 4: Hysteresis loop for the order parameter ci in a numerical simulation of 
the homogeneous kinetic equation with time varying p (see ( lA.ll) ). in dimension 2. 
The red curve is the theoretical curve, the blue one corresponds to the simulation. 
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Figure 5: Hysteresis loop for the order parameter ci in a numerical simulation of 
the homogeneous particle model with varying p (see (1A.2[) - (1A.3[) ). in dimension 2. 
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Figure 6: Free energy levels of the different equilibria (left), and difference of free 
energies of anisotropic and uniform equilibria (right), as functions of the density, in 
dimensions 2 and 3. The dashed lines in the right picture corresponds to unstable 
equilibria. At the density pi, the free energies of the stable anisotropic and the 
uniform equilibria are the same. 

distribution, while if p > pi, then the global minimum is reached for the family 
of stable von Mises-Fisher equilibria. However, there is no easy way to assess the 
value of pi numerically. We observe that the stable von Mises-Fisher distribution 
has larger free energy than the uniform distribution if p < pi and therefore consists 
of a metastable state. On the contrary, the uniform distribution has larger free en- 
ergy than the stable von Mises-Fisher distributions if p > pi and now, consists of a 
metastable state. 

The rates of convergence to the stable equilibria, following Proposition 13.41 and 
Theorem El are given by 

Ao = (n - 1)(1 ), for p < p^ = n, 

n 

A. = — — aJi- (--c-^^^^)j(1 + 2j)Y for p>p„ 

where Aq is the rate of convergence to the uniform distribution p, and A^ is the rate 
of convergence to the stable family of von Mises-Fisher distributions pM^Q, where k 
is the unique solution of the compatibility condition (I2.18|) such that n > n^,. Details 
for the numerical computation of the Poincare constant A^ are given in the appendix 
of [ID]. The computations in dimensions 2 and 3 are depicted in Figure[71 We observe 
that the rate of convergence to a given equilibrium is close to zero when p is close to 
the corresponding threshold for the stability of this equilibrium, and large when p is 
far from this threshold. Moreover, the rate A^ of convergence to a von Mises-Fisher 
distribution is unbounded as p — > oo, while the rate Aq of convergence to the uniform 
distribution is bounded by n — 1. 
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Figure 7: Rates of convergence to both types of stable equilibria, as functions of the 
density p, in dimensions 2 and 3. 

5 Macroscopic models, hyperbolicity, and difFu- 
sivity 

We now go back to the spatially inhomogeneous system. We want to investigate the 
hydrodynamic models that we can derive from the kinetic equation f l2.5p . 



5.1 Scalings 

In order to understand the roles of the various terms, it is useful to introduce di- 
mensionless quantities. We set to the time unit and xq = a to the space unit. 
We assume that the range of the interaction kernel K is R, meaning that we 
can write /^dxl) = -^K{'-^) (we recall that K is normalized to 1, so we still 
have /-ft'(|^|)d^ = 1). We also assume that K has second moment of order 1, 
i.e. K2 = 0{1), where 



K2 = -LKmm'd^- 



(5.1^ 



We now introduce dimensionless variables x = x/xq, t = t/tg, and we make the 
change of variables f{x,u,i) = Xq f{xoX,u,tQt), J'j{x,t) = XQj'f{xoX,tQt)/a. Fi- 
nally, we introduce the dimensionless quantities: 

V=^, H\Jj\) = H\Jf\)to, H\Jj\)=ri\Jf\)to. 

Xq 

In this new system of coordinates, the system (12.51) is written as follows (we still use 
the notation Jf{x,t) = J^ f{x,u,t) udu): 



dtf + io- VJ + z>(| JT/DV. ■ {P^^Qf f) = H\Jf\)A^f 
Jf{x,t) = {Kr!*Jf){x,t)dy, Uf = ^^, 



(5.2) 
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where we have dropped the tildes for the sake of clarity, and where K^, is the rescaling 
of K given by 

K,{x) = ^^K[^. (5.3) 

Now, by fixing the relations between the three dimensionless quantities fIS.ip . we 
define the regime we are interested in. We suppose that the diffusion and social 
forces are simultaneously large, while the range of the social interaction 77 tends 
to zero. More specifically, we let e -C 1 be a small parameter and we assume 
that f = 0{l/e) (large diffusion), v = 0{l/e) (large social force). In order to 
highlight these scaling assumptions, we define r", i^", which are all 0{1) and such 
that 

f = -t\ = -J. (5.4) 

Since r] is supposed to be small, using the fact that K is isotropic, we can first 
get the Taylor expansion of J7/ with respect to 77, using 05. 3p and 05. ip . when Jf is 
sufficiently smooth with respect to the space variable x: 

Jf = Jf + rflU^^Jf + 0(r/'). (5.5) 

Inserting this expansion into 05. 2p . and dropping all "hats" and "sharps", we are 
lead to: 

e{dtf + u ■ V.J) + K^rj' [v. ■ (P.x£;/) - m^A^/] = QU) + 0(r/^), (5.6) 

with 

QU) = -KI^/I)V. ■ (P^xf]; /) + ri\Jf\)A^f, 

■Jf 



Jf{x,t) = / f{x,u,t)udu, ilf 

Js ' K/ 1 

if = "^rrrPni ^.Jf + (% ■ a. J/) u'HJfDQf , 

where the primes denote derivatives with respect to \J\. We recover the same defini- 
tion of Q as in the spatially homogeneous setting 02. 8p , and the additional terms if 
and ruf do not depend on the velocity variable co (they only depend on Jf and its 
Laplacian A^Jf). 

Our plan is now to investigate the hydrodynamic limit e — )■ in this model, within 
two different regimes for the range of the social interaction 77: firstly, 77 = 0{e), and 
secondly, r] = 0{y/e). We have seen in [TT] that the second scaling allows us to 
retain some of the nonlocality of the social force in the macroscopic model, while 
the first one does not. Indeed, e corresponds to the characteristic distance needed 
by an individual to react to the social force, while 7/ is the typical distance at which 
agents are able to detect their congeners. The first scaling assumes that these two 
distances are of the same order of magnitude. The second one corresponds to a large 
detection region compared to the reaction distance. Which one of these two regimes 
is biologically relevant depends on the situation. For instance, we can imagine that 
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the first scaling will be more relevant in denser swarms because in such systems, far 
agents are concealed by closer ones. 

In both cases, we will write / as f^ to insist on the dependence on e. The limiting 
behavior of the function Z*^ as 5 — ?■ is supposed to be a local equilibrium for the 
operator Q, as can be seen in (15.61) . Keeping in mind the results of the previous 
section on the spatial homogeneous version, we will assume that /^ converges to 
a stable equilibrium of a given type in a given region. Depending on the type of 
equilibrium (uniform distribution or von Mises-Fisher distribution), we will observe 
different behaviors. 

5.2 Disordered region: diffusion model 

We consider a region where f^ converges as e — > to a uniform equilibrium p(x, t) 
which is stable. Therefore we must have p < pc- 

We first remark that we can integrate fl5.2p on the sphere to get the following 
conservation law (conservation of mass): 

dtPf^ + V, ■ Jf. = 0. (5.7) 

Therefore, if we suppose that the convergence is sufficiently strong, Jje converges 
to 0, and we get dtp = 0. 

To obtain more precise information, we are then looking at the next order in e in 
the Chapman-Enskog expansion method, in the same spirit as in the case of rarefied 
gas dynamics (see [S] for a review). We obtain exactly the same model as in [TU]. 
We prove the following theorem: 

Theorem 3. With both scalings r] = 0{^/e) and rj = 0{e), when e tends to zero, 
the (formal) first order approximation to the solution of the rescaled mean-field 
model 05.6J1 in a "disordered region" (where the solution locally converges to a stable 
uniform distribution) is given by 

r (x, CO, t) = p%x, t) - e ^^•V./>^(^^g (5.^ 



where the density p'^ satisfies the following diffusion equation 



[n — IjUTo \ 1 _ £- 



pc 



Proof We let p'^ = pfe and write f^ = p'^ + efl{x,uj,t) (so we have /§ /f dw = 0). 
The assumption is that /f is a 0{1) quantity as e — !■ 0. We then get 

Jf^=eJn. (5.10) 

and the model 05. 6p becomes: 

e{dtp' + u ■ V.p^) = -eu'{^)V^ ■ {P^^Jffp') + er^^Jl + 0{e^) + 0{rih). (5.11) 
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Additionally, using (15. 7p and f IS.lOp . we get that dtp"^ = 0{e). Therefore we can 
put edtp'^ into the 0{e'^) terms of (15.111) and get, in both scalings, 

A.A^ = -(V,p^ - u'{0){n - l)p'Jfi) ■ u + 0{e). 

To 

We can solve this equation for fl and, together with the fact that z^'(O) = ^^^ (thanks 
to the definition (I2.20p of pc) , we get 

/i = ( - ^^V.p^ + ^Jf.) ■ u + 0{e). 
^ ro(n-l) pc ^' 

This gives us, using (13.30 . that 

^ Ton{n - 1) Pc ' 

which implies that /f = nJfe ■ uj + 0{e) and that we have 

Jfi = 77^. —^^P' + ^(^) 



(n-l)nro(l-^ 

Therefore we obtain the expression (15. 8p of /^. Moreover, inserting this expression 
of J/e into the conservation of mass (15.70 gives the diffusion model (15. 9p . D 

Remark 5.1. The same remark was made in JT^ (see Remark 3.1 therein): the 
expression of ff, which is given by the 0{e) term of (15. 8p confirms that the approx- 
imation is only valid in the region where pc — p"^ ^ £■ The diffusion coefficient is 
only positive in the case where the uniform distribution is stable for the homogeneous 
model (p^ < pc) and it blows up as p^ tends to pc, showing that the Chapman- Enskog 
expansion loses its validity. 

5.3 Ordered region: hydrodynamic model and hyperbolicity 

We now turn to the derivation of a macroscopic model in a region where the local 
equilibria follow a given branch of stable von Mises-Fisher equilibria. More precisely, 
we suppose that the function /^ converges towards p{x,t)Mf^(^p(^x,t))Q,{x,t) in a given 
region, where n{p) is a branch of solutions of the compatibility equation (I2.22p 
defined for a given range of positive values of p, and which correspond to stable 
equilibria (in the sense of Theorem El). This implies that k is an increasing function 
of p. The goal is to prove the following theorem, which gives the evolution equations 
for p{x,t) and Q{x,t), assuming that the convergence of /^ is as smooth as needed. 

Theorem 4. We suppose that f^ converges as e — > towards p{x,t)MK(p(x,t))n{x,t), 
for a positive density p{x, t) and an orientation Q{x, t) G §, and where p H- K,{p) 
is a branch of solutions of the compatibility equation (12.220 . We also suppose that 
the convergence of f^ and of all its needed derivatives is sufficiently strong. Then p 
and VL satisfy the following system of partial differential equations: 

5tP + V, ■ (pcifi) = 0, (5.12) 

p {d,Vl + C2 (fi ■ V,)fi) + e Pn^V^p = /C2 5Pn^A,(pcifi), (5.13) 
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where IC2 is the scaling parameter corresponding to the limit of K2— as e ^ 0, and 
where the coefficients ci, C2, Q and S are the following functions of p{x,t) (where 
the dependence on p for k or on K,{p) for c, c or j is omitted when no confusion is 
possible): 

Clip) = c(«:) (5.14) 

(cos^/i«(cos^)sin2e)M _ Jq cos9 h^{cos9) e''^°'^ sin^^d^ 
C2(p) - c{K) - (^h^^cose)sm^e)M ~ J^h^{cose)e---^<^ sin'^ede ' ^ ' 

Q{p) = - + {c-c)^— = —, '- ^^, 5.16 

Kj2 /^- 1 

C ^ K 

The function h^ is defined below at Proposition 15.51 



.W = M(!i^ + g). (5.17) 



Proof. The first equation fl5.12p (for tlie time evolution of p) can easily be derived 
from the conservation of mass (15.71) . since the hypotheses imply that Jfs converges 
to pc{K,{p))fi, and thanks to (15.51) . we obtain: 

9tP + V,-(pc(K(p))fi) = 0. 

The main difficulty is the derivation of an equation of evolution for Q, since we do 
not have any conservation relation related to this quantity. To this aim, the main 
tool consists in the determination of the so-called generalized coUisional invariants, 
introduced by Degond and Motsch [12] to study the model corresponding to the 
case u = T = 1 and the scaling 77 = 0{e) in our setting. These generalized coUi- 
sional invariants were then used successfully to derive the same kind of evolution 
equation for some variants of the model we are studying (see [lU] when z^(| J|) = | J| 
and T = 1, [IS| for the case where 1/ is a function of p, and where the interaction 
is anisotropic, [13] for another type of alignment, based on the curvature of the 
trajectories, and [TT] in the case of the second scaling rj = O(y^) when u = t = 1). 
The idea is to introduce, for a given k > and ^2 G §, the operator L^q (linearized 
operator of Q) : 



L^nif) = A J - /.V^ ■ (/P^xfi) = V, 



M^nvJ ^ 



M. 



kQ , 



SO that we have Q{f) = r(| Jj|)Lfc(|j^|)Q^(/). And we define the set of generalized 
coUisional invariants C^q: 

Definition 5.1. The set C^n of generalized coUisional invariants associated to k, &M. 
and fi G S, is the following vector space: 

C^n = i ^1 / L^nif) tpduj = 0, V/ such that Pn±Jf = o\ . 

Hence, if ■?/; is a coUisional invariant in C^q, we have J^^gQ{f)ipdu} = for any 
function / such that k{\Jf\) = k and flf = Q. 

The determination of C^n has been done in [TH]. We recall the result here: 
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Proposition 5.1. There exists a positive function h,^ : [—1, 1] — !■ M such that 

C^^ = {uj^ h^{Lo ■ n)A ■ w + C, C e M, A e M", with A-n = 0}. 

More precisely, /Ik(cos6') = ^^^, where (?«; is the unique solution in the space V 
of the elliptic problem L*^g{9) = sin 6', where 

Llg{e) = -sin2-"ee-''=°'^^^(sin"-2^e'''=°^V(^)) + ^9i^), 

We now define the vector-valued generalized collisional invariant associated to k 
and Q as 

^kq(w) = K{u ■ Q)Pq±U, 

and we have the following useful property: 

V/ such that k{\Jf\) = k and Vlf = Vl, / Q{f) ipno. dw = 0. 

The next step consists in multiplying the rescaled kinetic model (15.61) by ^V^K^f^/e, 
with K^ = k{\Jfe\), and to integrate it on the sphere. We get, 

2 4 



P(f,.)x (X^ + K^^y + Z']) = 0[-L 



where 



2 

Now we can pass to the limit £ — >• 0. We denote by /C2 the limit of -ft'2^; which 
makes sense in both scalings (either tj = 0{e), and /C2 = 0, or 77 = 0{^/e) and we 
suppose that /C2 is a positive quantity), and we obtain 

Pn^ {X + }C2[Y + Z]) = 0, (5.18) 

where, since we suppose that /^ — ;■ p(x, t)MK(p(a;,t))n{a;,t)), we have 

X= f (dtipM^n) + CO ■ V,{pM^n)) K{uj ■ Q) u du , 

Y= f V^ ■ iP^±ipM^^pM^n) K{oo ■n)ujduj, 

Z = I mpM^a^u{pMf,n) hf,{uj ■ (]) w dw . 

The computation of Pq± X has been done in [IS]: we get 

Pn^X ={h^{cos 9) sin^ 9) MP^^ (a^fi + 5(f] ■ V.)^) + 6 Ff,xV,.p, (5.19) 

n — I 
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where c and G are given by f l5.15p and the first expression of fl5.16p . We now 
compute Pq±Y and Pfi±Z. Since V^j ■ {P^±A) = —{n — l)A-uj for any vector A e M", 
we get 

Writing u = cos9Q + sin6'f (orthogonal decomposition with v G Ein-2), and using 
the fact that /s^,, f dw = and /§^ ^ 'i^ 8) v dv = ^^^y-^f^-^ ' ^^ obtain 

Pn^ Y = - {h^{cos 0) sin^ 9)m P — Pn^ipM^n 

n — 1 

= - {K{cos 9) W 9) Mv{pc) ^ - 1 +^^ p^^A^(pc(]). (5.20) 

[n — l)c 

Finally, since A^(MKn) is a function of cos^, the same decomposition and argument 
shows that we have Pq,^Z = 0. Inserting (15.191) and (15.201) into (15.180 and dividing 
by ^^(/ik(cos6') sin^ 9)m ends the derivation of (15.131) . with 6 given by (I5.17p . 

We finally derive the expression of 6 given by the right-hand side of (15.160 . 
We differentiate the compatibility condition pc = j with respect to k, (in a given 
local branch of solutions), and we get c^ + p^ — '^- ^s was shown in (TU], we 
have ^ = 1 — (n — 1)^ — c^, therefore we get 

K dp K dp K dj K,^ , ,C 9, , ^ Kdj K , 

- / = c - / = -/ - p -(1 - (n - 1)- - c^) = (n - 1 + -/ - - + Kc), 
p dn J dn J dn J k j dn c 

and finally, thanks to the first expression of (15.160 . we have 

1 c-c n-^ + KC-l + ^^ 

n n-l + fg-f + «:c ^(n - f + kc - 1 + f £) ' 

which gives an expression of in terms of k and the functions c, c and j only. D 

5.4 Hyperbolicity 

As shown in [TT |[T^[TB] . and more precisely in ^U\ when the coefficients ci, C2 and G 
depend on p, we have the following 

Proposition 5.2. The self- organized hydrodynamic (SOH) model (I5.12p -( 15.13p is 
hyperbolic if and only if Q > 0. 

In that case it has been proved in [TT] that the SOH model is locally well-posed 
in dimension 2 (provided 5^0) and in dimension 3 (for the scaling where /C2 = 0, 
with an additional condition for the orientation of the initial data). Therefore, in 
this section, we study the sign of these coefficients in some generic situations. 

Conjecture 5.1. For all k > 0, we have < c{k) < c{k). Consequently, in the 
SOH model, we have 6 > 0, and the SOH model is well-posed if it is hyperbolic. 



44 



Numerically, this conjecture is clear, at least in dimension 2 and 3, as can be 
seen in Figure 3 of [lU]. We know it is true when k is small or large, thanks to the 
asymptotics of c and c given in [15] : 



l^ - n^)^"" + 0('^') as «: -> 0, 



(5.22) 



^~\. n+1 I (n+l)(3n- 7) , ^/ -Sn l^O.ZOJ 

2k "•" 24k2 



These asymptotics can also help us to know if the system is hyperbolic in various 
regimes. In the next four propositions, we provide different cases where we can deter- 
mine the hyperbolicity of the SOH model with simple assumptions on the behavior 
of the function k. The first result is about non-hyperbolicity in the neighborhood of 
the critical threshold p^, for a first order phase transition. 

Proposition 5.3. Suppose that there is a first order phase transition with hystere- 
sis as described by Proposition \4-^ If Conjecture 15. Jl is true, then the SOH model 
associated to the branch of stable von Mises-Fisher equilibria (for p > p^) satis- 
fies 6(p) < if p is sufficiently close to p*. The SOH model is not hyperbolic. 

Proof We have (^)'(k*) = and (^)'(fi:) > for k > k^. This gives that ^ —^ +oo 
as K — !■ K*, and then we use ( I5.16P and Conjecture lS.ll to get that 6 — )■ — oo as k — )■ k^, 
(for n > /€*). D 

We now provide the same type of proposition in the neighborhood of the critical 
threshold pc in the case of a second order phase transition. The following proposition 
gives a strong link between hyperbolicity and the critical exponent of a second order 
phase transition: it is hyperbolic when the critical exponent (3 is greater than |, and 
not hyperbolic when /3 < | (this threshold value | also corresponds to the lowest 
possible critical exponent which can appear in the case of the simple criterion given 
by Lemma [3]). 

Proposition 5.4. We suppose that there is a second order phase transition as de- 
scribed by Proposition \4-l\ and we consider the SOH model associated to the von 
Mises-Fisher equilibria (for p > p^). We suppose furthermore that there is a critical 
exponent (3 as stated in Definition \4.1\ and we assume that this estimation is also 
true at the level of the derivative: 

-— ~ naof3(p - Pc)^'-^, as p^ pc- 
dp 

Then 

(i) If P < ^, then Q{p) < if p is sufficiently close to p^. The SOH model is not 
hyperbolic. 

(a) If (3 > \, then 6(p) > if p is sufficiently close to pc. The SOH model is 
hyperbolic. 
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(Hi) If (3 = ^ and Oq 7^ \ 5""*" > then when p is sufficiently close to pc, ©(p) is of 
the sign of ^^^^ - a^. 

Proof. We have k{p) ~ nao{p — PcY , as p -^ pc- So we get ^^ ~ /3-r^- Finally, 
using fl5.22p - fl5.23p . we get c — c ^^ ~ 2n(n+2) '^ as k — )► 0. We can then obtain an 
equivalent of 6 as p —^ Pc, with fl5.16p : 

(^ ^PcO'o] 1 if /5 = 1 

where the last expression is valid only if ^^ 7^ ifn+2) ■ "^^^ ^^S^ '^^ ® ^^ then directly 
given by these equivalents, and this ends the proof. D 

It is possible to refine Proposition 14.21 in order to have the critical exponent 
estimation on the level of the derivative, and then express the hyperbolicity of the 
system with the help of the expansion of k only. In summary, we get the following 
proposition, the proof of which is left to the reader: 

Proposition 5.5. If k satisfies: 

71 

k'(\J\) = a(q + 1)1 Jl" + o(\.m as \J\ -^ 00, 

Pc 

then we have 

(i) if q < 2 and a > 0, the critical exponent is given by (3 = - and the SOH model 
is hyperbolic when p is sufficiently close to pc- 

(a) if q = 2 and a > 4 3fa-L2) > then P = ^ and the SOH model is hyperbolic when p 
is sufficiently close to pc- 

(Hi) if q = 2 and — o " , „, < a < , J} , „^ , then = k OLnd the SOH model is not 
hyperbolic for p close to pc. 

We finally give a result about hyperbolicity when p is large, depending on the 
behavior of A; as | J| — > 00. 

Proposition 5.6. We suppose that k{\J\) ~ a \J\^ as \J\ -^ 00 (with a,b> 0), and 
that this equivalent is also true at the level of the derivative: k'{\J\) r^ ab \J\^~^ . We 
consider the SOH model associated to a branch of stable von Mises-Fisher equilibria. 

(i) If < b < 1, then for p sufficiently large, Q{p) < and the SOH model is not 
hyperbolic. 

(a) If b > 1, then for p sufficiently large, 0(p) > and the SOH model is hyper- 
bolic. 
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(in) Ifb= 1, we have to make stronger hypotheses on k. For example, if we suppose 
that, as I J| — ;■ oo, we have k{\J\) = a\J\ + r + o(l) and k'{\J\) = a + o(| J|~^) 
with r 7^ ^^, then for p sufficiently large, Q{p) is of the sign of r — ^^. 

Proof. We first use the expansion (15.221) to get that n— 1 — ^ + fi;c~ —^^ as k — ;■ oo, 
and that c — c = —^ + ^^ + o(fi;~^). Using Hypothesis 12.21 the assumptions 
become j{k) ~ (f)i and £ = [k'{j{K))]~^ ~ {aby^^ft"^ as k ^ oo. This 
gives p = 4^ ~ (^)6 as K — ^ oo, which can be inverted to get k ~ ap^ as p —^ oo. 
Finally, for b y^ 1, we get, with the left part of (I5.2ip : 

e(p) ~ (l - -) — - as p ^ +00. 
^ bap" 

This proves the first two points. In the case where 6 = 1, we suppose that we have 
the expansions A;(| J|) = a\J\ + r + o(l) and k'{\J\) = a + o{\J\~^) as \J\ — )■ oo. 
Then j{k) = ^{k — r) + o(l) and ^f = ^ + o{k~^). And we finally get, using the left 
part of (15.211) . as k -> +oo: 



e = - 



1 -K-' + ^K-^ + 0{k-^) 



Since k ~ ap as p —^ oo, we have 6(p) ~ ^(r — ^^)p~'^ as p — ^ oo and this proves 
point (iii). D 

Remark 5.2. The case 6 = can also be treated if we assume Conjecture \5.1[ 
This corresponds to the case where k takes values on [0,Kmax) with i^max < oo. // 
furthermore we assume that its derivative satisfies k' ~ a| J|~^ (with b > 1 and a > 0) 
as \J\ -^ oo, then after the same kind of computations we get that 0(p) — )■ +oo 
as p ^ +00, and the system is hyperbolic. 

Let us now comment these results in the case of specific examples. The case 
where z^(| J|) = \J\ and t(| J|) = 1 corresponds to the model studied in [TU]. It was 
shown to be non hyperbolic (numerically for all p > pc = n, and theoretically for p 
large or close to pc). We now see that it corresponds to points (iii) of Proposition EIS 
and Proposition 15.51 which are the special cases separating hyperbolicity to non- 
hyperbolicity. A really slight change in the function k in this model could easily 
lead to hyperbolicity, while nearly keeping the same phase transition phenomena, 
from the point of view of equilibria. 

The case studied in Section 14.31 and leading to a first order phase transition 
corresponds to the function /i;(|J|) = |J| + |Jp. Thanks to Proposition 15.31 and 
to 15.61 we get that the corresponding SOH model is not hyperbolic in both regimes: 
when p is close to p* and when p is sufficiently large. Numerical computations of 
the coefficient suggest that this is the case for all the values of p > p* (at least in 
dimensions 2 and 3). 

Finally, we are interested in the original model presented in [12], where u and r 
are constant. We remark that hypotheses 12.11 and 12.21 do not cover this model, 
but we can see it as a limiting case of a regularized u satisfying such hypotheses, 
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such as z^^(|J|) 
Proposition 



= J,<j< - In that case, we have pi = nero, and by Lemma E] and 
we get that there is a second order phase transition with critical 
exponent 1. Furthermore, with Remark 15.21 and Proposition 15.41 we get that the 
corresponding SOH model is hyperbolic when p is large or close to pc- Figures [S] 
and El correspond to the plots of the phase diagram (the order parameter c) and of 
the function G for three different values of e, with tq = | and n = 2. We observe 
that the system is always hyperbolic. 
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Figure 8: Order parameter ci, as function of the density p, in dimension 2, for the 
regularized model. 



We get the same conclusion for a regularization given by z^^(|t/|) 
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with a critical exponent f3 
sit ion (15.51) if Tq > 



= I this time, and k satisfies the condition (ii) of Propo- 
. =. This gives a practical example of a second order phase 

y/2n(n+2) 

transition with the minimal critical exponent such that the associated SOH model 
is hyperbolic in the neighborhood of the threshold pc (indeed, in that case, thanks 
to Proposition (15.41) . we must have /? ^ |). 



6 Conclusion 

In this work, we have provided a comprehensive and rigorous description of phase 
transitions for kinetic models describing self-propelled particles interacting through 
alignment. We have highlighted how their behavior results from the competition 
between alignment and noise. We have considered a general framework, where both 
the alignment frequency and noise intensity depend on a measure of the local align- 
ment. We have shown that, in the spatially homogeneous case, the phase transition 
features (number and nature of equilibria, stability, convergence rate, phase diagram. 
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Figure 9: Coefficient 9, as function of the density p, in dimension 2, for the regu- 
larized modeL 

hysteresis) are totally encoded in the function obtained by taking the quotient of the 
alignment and noise intensities as functions of the local alignment. The phase transi- 
tions dealt with in this paper belong to the class of spontaneous symmetry-breaking 
phase transitions that also appear in many physics systems such as ferromagnetism, 
liquid crystals, polymers, etc. We have also provided the derivation of the macro- 
scopic models (of hydrodynamic or diffusion types) that can be obtained from the 
knowledge of the stable equilibria and classified their hyperbolicity. In particular, 
we have provided a strong link between the critical exponent in the second order 
phase transition and the hyperbolicity of the hydrodynamic model. In the future, 
we will investigate how the hydrodynamic and diffusion regimes can be spatially 
connected through domain walls and find the dynamic of these domain walls. 

A Appendix: numerical methodology for the hys- 
teresis simulation 

In this appendix, we give more details on the computation of the hysteresis loop 
provided in section 14. 3[ In order to highlight the role of the density p as the key 
parameter for the phase transition, we introduce the probability measure / = - and 

we rewrite the homogeneous kinetic equation 02.71) in terms of /. We get 



dtf = Tip\Jj\)Aj- z/(p| Jr|)V^ ■ iP^.njf). 



(A.l) 



When p is constant, this equation is equivalent to (12.71) . We will now consider p 
as a parameter of the equation (and not anymore a parameter for the mass of 
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initial condition, since / is now a probability measure), but the long time behavior 
(equilibria, stability, convergence) is still given by this parameter p. Finally, we let p 
vary slowly with time (compared to the time scale of convergence to equilibrium, 
given by Figure [7]), as we expect it would be the case in the spatial inhomogeneous 
framework given by the kinetic equation (12. 5p . 

A.l Simulation at the kinetic level. 

Let us now present how the numerical simulations of the system flA.ll) in dimen- 
sion n = 2, depicted in Figure HI have been obtained. We start with an initial 
condition which is a small perturbation of the uniform distribution, and we take a 
varying parameter of the form p = 1.75 — 0.75cos(^t), with T = 500. We use a 
standard central finite different scheme (with 100 discretization points), implicit in 
time (with a time step of 0.01). The only problem with this approach is that the 
solution converges strongly to the uniform distribution for p < p^. So after pass- 
ing pc, the linear rate of increase for Jj is given by — — 1, by virtue of (13. 9j) . and 
is very slow when p is close to pc- So since Jj is initially very small when passing 
the threshold p = pc, the convergence to the stable von Mises-Fisher distribution 
is very slow. Two ideas can be used to overcome this problem: either injecting 
noise in the system, or more efficiently, adding a threshold e and strengthening | Jj\ 

when 11/ — l||oo ^ £, replacing / at the end of such a step by 

f + max{0,e- \\f - l\\^)nj-uj. 

We note that after this transformation, we still have ||/ — l||oo ^ e if it was the case 
before applying the transformation. 

Figure m depicts the result of a numerical simulation with a threshold e = 0.02. 
We clearly see this hysteresis cycle, which agrees very well with the theoretical 
diagram. The jumps at p = p,,. and p = pc are closer to the theoretical jumps 
when T is very large. 

A. 2 Simulations at the particle level. 

Now, since the kinetic equation (12.51) comes from a limit of a particle system, we are 
interested in observing this hysteresis phenomenon numerically at the level of the 
particle system, where noise is already present in the model, since it is a system of 
stochastic differential equations. 

As for (12. 5p . it is easy to derive the mean-field equation (lA.ip . in the spirit of 
Proposition 12.11 from the following system: 



dUi = (/(pI J|)P„_i SJd( + v'2r(p| J|)P„_i o cIB;, (A,2) 

Here, once again, the parameter p is a parameter of the equation, which can be 
variable in time. We perform numerical simulations of this system for a large number 
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of particles, with p varying as in the numerical simulation of the kinetic model. As 
before, we start with a initial condition which consists of A^ = 10000 particles 
uniformly distributed on Si, and we take p = 1.75 — 0.75 cos(|:t), with T = 500. We 
use a splitting method for the random and the deterministic parts of this equation 
(with a time step of 0.01). We then plot the order parameter c, given by \J\. The 
result is given in Figure O 

Let us remark that, thanks to the central limit theorem, the mean J oi N vectors 
uniformly distributed on the circle has a law equivalent to a multivariate normal 
distribution in M^ centered at 0, and with covariance matrix j^l2- Therefore \J\ is 
equivalent to a Rayleigh distribution of parameter -7=, and so the mean of \J\ is 

equivalent to ^^- In our case, that gives a mean of \J\ of approximately 0.009, of 

the same order as in the previous section, since the threshold e ensures that, when / 
is close to the uniform distribution | J~| Ri | with e = 0.02. 
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